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ultrasonic  mass  fuel  flowmeter.  The  report  consists  of  a  discussion 
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INTRODUCTION 


From  both  diagnostic  and  control  viewpoints,  it  is  important  to  know  ac¬ 
curately  the  value  of  the  fuel  mass  flow  rate  Nt.  Since  the  heating  values 
for  typical  fuels  such  as  JP-4,  JP-5,  diesel  fuel  and  Avgas  grades  from 
80  to  145  octane  are  within  about  1%  of  1 9 ,  000  Btu/ib,  it  is  clear  that  the 
maximum  power  that  can  be  extracted  from  the  fuel  is  proportional  to 

In  engine  diagnostic  studies,  one  is  interested  in  measuring  engine  per¬ 
formance  from  start-up,  at  relatively  low  flow  rates  of  about  50  1b/hr,  up 
to  maximum  running  conditions,  where  Id  =  2000  lb/ hr  in  some  engines. 
Thus,  a  range  of  ~ 40 : i  is  of  interest.  For  diagnostics,  response  time 
of  500  ms  presently  appears  to  be  adequate.  Accuracy  of  better  than  1% 
is  desirable  at  high  flow  and  ~5%  at  the  low  end. 

In  engine  control  systems,  one  can  control  in  terms  of  the  temperature  of 
the  combusted  fuel  and/ or  in  terms  of  M.  The  same  Id  ranges  and  accu¬ 
racy  apply  as  for  diagnostics,  but  response  time  should  be  in  the  10  to  30 
ms  range. 

While  a  number  of  commercially  available  flowmeters  could  meet  most 
or  all  of  the  above  requirements,  there  are  additional  requirements  im¬ 
posed  by  the  situations  where  Army  aircraft  gas  turbine  engines  art; 
tested  and/  or  used.  For  example,  the  fuel  composition  maybe  an  un¬ 
known  mixture  of  two  or  more  fuels.  The  fuel  temperature  may  range 
between  wide  limits.  The  fuel  may  be  contaminated  by  different  types  of 
particles  of  various  sizes,  and  by  corrosive  material.  The  engine 
environment  includes  noise  and  vibration.  The  flowmeter  measuring 
cell,  transducers  and  electronics  should  be  small  and  lightweight,  and 
rugged  enough  to  withstand  the  expected  operating  conditions. 

The  main  problem  with  previous  (turbine)  flowmeters  has  been  clogging 
within  ten  hours  under  contaminated  fuel  flow.  The  present  ultrasonic 
flowmeter  employs  recessed,  nonintrusive  transducers  and  has  no  ro¬ 
tating  parts.  Accordingly,  it  is  not  subject  to  clogging  and  has  continued 
to  operate  after  being  contaminant-tested  for  approximately  one  hundred 
hours,  i.  e,  ,  substantially  longer  than  specified  in  MIL/-F-861  5. 

A  logical  development  of  the  present  breadboard  system  is  expected  to 
meet  accuracy  and  range  requirements.  Contractual  objectives  for  upper 
flow  range  and  speed  of  response  have  already  been  exceeded  by  factors 
of  2.  5  and  25,  respectively. 
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The  following  sections  of  this  report  define  the  problem  of  measuring 
mass  flow  rate  in  greater  detail,  indicate  the  background,  explain  the 
ultrasonic  approaches  to  the  problem,  and  describe  the  measuring  and 
testing  methods  and  results,  Conclusions  are  based  on  our  analysis, 
measurements,  tests,  and  results;  and  recommendations  are  offered 
relative  to  achieving  the  remaining  objectives. 
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STATEMENT  OF  THE  PROBLEM 


The  purpose  of  this  program  was  to  develop  an  ultrasonic  mass  flowmeter 
to  meet  the  following  objectives: 

a.  Operate  over  the  temperature  range  of  -65°F  to  +160°F 
( -  54°C  to  +7 1  °C). 

b.  Measure  fuel  flcv<  over  the  range  of  50  lb/hr  to  2000  lb/ hr, 
with  an  accuracy  of  2.  5  lb/ hr  at  a  flow  of  50  lb/ hr,  and 

10  lb/hr  at  a  flow  of  2000  lb/hr,  with  linear  variation  be¬ 
tween  this  range. 

c.  Maximum  integration  time  of  0.  5  second. 

d.  The  transducer  was  to  have  no  moving  parts. 

e.  The  transducer  was  not  to  restrict  or  obstruct  the  flow. 

f.  The  transducer  was  to  be  sized  and  configured  so  as  to  be 
capable  of  being  mounted  on  an  Army  aircraft  gas  turbine 
engine. 

In  designing  an  ultrasonic  flowmeter,  one  must  translate  these  objectives 
into  ultrasonic  terms.  P'or  example,  items  (a)  and  (f)  translate  to  trans¬ 
ducer  materials  selection.  Barium  titanate,  in  common  use  below  its 
Curie  point  of  about  120  C,  would  not  necessarily  be  appropriate  because 
of  the  higher  ambient  temperatures  encountered  at  no-flow,  after  engine 
turnoff.  Of  the  many  other  transducer  materials  available,  the  final 
choice  may  depend  on  variation  of  transduction  and  impedance  values 
versus  temperature.  As  a  further  example,  item  (h)  translates  to  ultra¬ 
sonic  measurement  of  flow  velocity,  sound  speed  in  the  fuel,  and  density, 
with  calculable  error  limits. 


BACKGROUND 


The  deficiencies  of  the  mass  flow  rate  sensors  currently  used  to  diagnose 
Army  aircraft  gas  turbine  engine  performance  include: 

•  Turbine-type  sensors  erode,  deg.  '.e,  cease  to  operate  due  to 
dirt  -  lifetime  limited. 

•  High  accuracy  and  high  reliability  not  found  in  same  sensor. 

•  Sensitive  to  vibration. 

•  Calibration  depends  on  fuel  type/ composition. 

•  Density  correction  depends  on  factors  besides  density. 

•  Rotating  turbine-type  sensor  must  be  maintained  inside  the  fuel 
pipeline.  Repair  requires  disassembly,  downtime. 

Recognizing  the  significance  of  the  above  problems,  and  the  need  for  an 
approach  sufficiently  new  and  advanced  so  that  it  could,  by  design,  sub¬ 
stantially  avoid  these  problems,  a  program  was  initiated  in  June  1971  to 
develop  an  ultrasonic  mass  flowmeter  system. 

Historically,  the  fundamental  acoustic  principles  underlying  most  present 
ultrasonic  flowmeters  are  traceable  to  Doppler  (the  "Doppler  effect")  or 
to  Newton,  who  recognized  that  the  time  for  sound  waves  in  air  to  propa¬ 
gate  between  two  points  depended  on  the  wind  velocity.  Worldwide,  over 
one  hundred  ultrasonic  flowmeters  presently  measure  liquid  flow,  and 
occasionally  gas  or  particulate  flow,  in  industrial  installations.  Additional 
ultrasonic  flowmeters  are  used  in  biomedical  applications,  for  measuring 
blood  flow  or  the  motion  of  organs.  The  Doppler  effect  is  utilized  in  most 
biomedical  flowmeters,  while  the  transit  time  difference  (upstream  minus 
downstream)  is  utilized  in  most  industrial  flowmeters. 

In  this  program  we  pursued  initially  the  Doppler  effect,  or  scattering  ap¬ 
proach,  to  measure  flow  velocity,  but  later  changed  to  a  transit  time,  or 
phase  shift,  approach.  Reflection  principles  were  also  utilized,  to  mea¬ 
sure  fuel  density. 
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ULTRASONIC  APPROACHES  TO  THE  PROBLEM 


As  stated  on  page  2,  part  of  the  problem  translates  to  ultrasonic  mea¬ 
surements  of  fuel  flow  velocity,  sound  speed  in  the  fluid,  and  fuel  density. 
Table  I  lists  the  specific  ultrasonic  approaches  that  were  used  for  these 
three  measurements. 


TABLE  I. 

ULTRASONIC  APPROACHES  USED  IN 

MEASURING  FUEL  FLOW  VELOCITY, 

SOUND  SPEED,  AND  DENSITY 

Fuel  Parameter 

Measurement  Approach 

Flow  Velocity 

Doppler  scattering;  phase  difference  on 
transmission,  using  phase  meter, 

Sound  Speed 

Transit  time  of  transmitted  pulse,  using 
intervalometer , 

Density 

Reflection  coefficient  at  probe/ fuel  inter¬ 
face,  using  special  probe  -nd  circuitry. 

To  clarify  these  approaches,  we  explain  briefly  the  basic  ideas  underlying 
acoustic  scattering,  transmission  and  reflection  measurements.  Then  we 
present  a  more  detailed  description  of  the  three  tasks  into  which  the  pro¬ 
gram  was  divided:  Task  I,  Theoretical  and  Analytic  Optimization;  Task 
II,  Design  and  Fabrication;  and  Task  III,  Tests. 

SCATTERING 


Doppler  flowmeters  measure  the  change  in  frequency  of  a  wave  scattered 
off  moving  particles  in  the  fluid.  The  Doppler  shift  is  proportional  to  the 
ratio  of  particle  velocity  to  sound  speed.  If  the  particles  are  small 
enough  each  acquires  the  local  fluid  velocity.  If  they  are  uniformly 
distributed  and  uniformly  insonified,  the  scattered  wave  can  be  processed 
to  yield  the  average  flow  velocity.  When  the  scatterers  are  not  uniformly 
distributed,  range-gated  Doppler  techniques  are  sometimes  employed  to 
obtain  the  flow  profile,  from  which  the  average  flow  velocity  may  be 
computed. 

The  "particles"  may  be  foreign  bodies  such  as  contaminants,  c  r  bubbles, 
or  they  may  be  regions  of  fuel  having  an  acoustic  impedance  slightly  dif¬ 
ferent  from  the  neighboring  fuel.  In  Doppler  blood  flowmeters,  the  red 
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corpuscles  are  the  particles  which  scatter  ultrasound  having  a  frequency 
of  about  5  to  10  MHz.  However,  Doppler  scattering  has  also  been  ob¬ 
served  in  pure  distilled  liquid  metals  such  as  sodium-potassium  (NaK). 

1  he  Doppler  approach  was  selected  at  the  beginning  of  this  program, 
because  at  that  time  it  appeared  to  offer,  at  least  theoretically,  the  best 
way  of  measuring  average  flow  velocity  with  minimum  disturbance  of 
simple  pipelines.  Calculations  of  the  power  scattered  at  20  MHz  (wave¬ 
length  approximately  0.  05  mm,  or  50  pm)  from  fuel  contaminated  accord¬ 
ing  to  MI.L-E-5007C  showed  that  scattering  should  be  adequate  even  in 
fuels  containing  approximately  1%  of  the  specified  contaminant.  However, 
tests  on  uncontaminated  liquid  did  not  provide  enough  scatter  to  be  de¬ 
tected  with  our  equipment. 

The  theory  of  scattering  of  sound  in  a  fluid  by  obstacles  which  are  rigid 
or  nonrigid  (including  bubbles)  has  been  well  developed  for  simple  shapes 
such  as  spheres  or  cylinders.  Without  tracing  the  theory  all  the  way  back 
to  Rayleigh,  we  merely  present  one  of  his  simple  but  important  theoret¬ 
ical  results,  for  a  spherical  scatterer  whose  radius  is  very  small  com¬ 
pared  to  wavelength.  This  result  is  that,  for  so-called  Rayleigh  scattering, 
the  scattered  power  is  directly  proportional  to  the  sixth  power  of  the  radius, 
and  inversely  proportional  to  the  fourth  power  of  the  wavelength.^  -  2) 

This  Rayleigh  scattering  result  suggests  that  one  could  increase  the  scat¬ 
tered  power  by  increasing  the  frequency  (decreasing  the  wavelength). 
However,  the  sound  attenuation  coefficient  increases  as  the  square  of 
frequency.  As  a  practical  compromise  on  frequency,  we  calculated  that 
for  a  path  of  5  cm,  the  maximum  frequency  would  be  30  MHz,  with  10  to 
20  MHz  probably  being  the  most  appropriate  range  to  utilize  (Appendix  I). 

We  also  considered,  but  did  not  pursue,  the  possibilities  of  increasing 
acoustic  scattering  by  increasing  the  size  (and  number)  of  scatterers, 
either  by  introducing  air  bubbles  or  by  cavitating  the  fuel. 

At  this  point  in  the  program,  a  new  theoretical  development  unfolded, 
which  predicted  that  the  average  flow  velocity  could  be  determined  by  a 
new  transmission  measurement.  This  is  discussed  next. 

TRANSMISSION 


Transit  time  or  phase  shift  flowmeters  usually  measure  the  difference  in 
transmission  time  or  phase,  with  oblique  propagation  upstream  and  down¬ 
stream^  These  observed  differences  are  proportional  to  the  average  flow 
velocity  vft  in  the  acoustic  path,  divided  by  the  square  of  the  sound  speed, 
C2.  Previous  techniques  to  account  for  nonuniform  flow  velocity  profile 


include  iterative  methods  based  on  estimates  or  knowledge  of  the  Reynolds 
number,  and  Gaussian  quadrature  methods  which  compute  the  average 
flow  from  measurements  across  four  parallel  chords.  Neither  of  these 
earlier  approaches,  however,  sampled  all  of  the  fluid. 

A  novel  way  of  obtaining  the  average  flow  velocity  was  utilized  in  the 
present  program.  An  ultrasonic  measuring  cell  was  designed  such  that 
all  the  flowing  liquid  could  be  interrogated.  The  interrogation  was  ac¬ 
complished  by  a  pair  of  5  MHz  continuous  plane  waves  which  ideally  re¬ 
main  substantially  undistorted  in  the  flowing  liquid.  To  separate  upstream 
from  downstream  waves,  each  wave  was  phase-coded  witha  pseudo-random- 
noise  code. 

Since  the  observed  differences  in  transmission  time  are  proportional  to 
c  +  v,  where  c  >>  v,  it  is  clear  that  a  small  fractional  error  in  c  is 
equivalent  to  a  large  fractional  error  in  v.  To  minimize  this  error  con¬ 
tribution,  the  upstream  and  downstream  paths  should  be  as  identical  as 
possible.  In  practice,  therefore,  it  is  desirable  to  use  the  same  pair  of 
transducers  for  both  upstream  and  downstream  measurements.  This  was 
done  in  the  present  design. 

To  eliminate  the  dependence  of  the  output  signal  on  c^,  one  can  arrange  to 
measure  the  difference  between  two  frequencies,  each  frequency  being  in¬ 
versely  proportional  to  one  of  the  travel  times.  This  was  not  done  in  the 
present  program,  because  at  the  time  the  Doppler  approach  was  abandoned, 
the  two-variable- f requency  technique  required  to  continuously  measure  v 
independent  of  c  had  not  yet  been  conceived.  In  any  event,  the  present 
v/  c^  measurement  would  still  be  a  logical  first  step. 

In  passing  we  may  note  that,  compared  to  the  transmission  approach,  a 
potential  advantage  of  the  Doppler  approach  to  mass  flow  metering  is  that 
ihe  Doppler  shift  is  proportional  to  v/c.  Multiplying  this  quotient  by  the 
pc  output  of  the  reflection  coefficient  densitometer  yields  an  output  pro¬ 
portional  to  W,  without  the  need  to  measure  c  explicitly.  We  may  further 
note  that,  alternatively,  if  one  can  use  a  densitometer  approach  that  di¬ 
rectly  yields  p  alone  (such  as  certain  resonant  structures  whose  resonant 
frequency  is  a  function  of  p  ,  independent  of  c)  then  a  transmission  mea¬ 
surement  of  upstream  and  downstream  frequencies,  whose  difference  is 
proportional  to  v,  appears  most  appropriate.  Multiplying  these  latter  out¬ 
puts  fields  M  =  p  v. 

In  this  program  a  separate  transmission  path  normal  to  the  fuel  flow  di¬ 
rection  was  utilized  to  measure  a  time  interval  T,  which  is  inversely 
proportional  to  c,  and  independent  of  v,  with  a  precision  of  _+  0.  05%. 
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REFLECTION 

The  reflection  approach  to  measuring  the  properties  of  a  medium  is  based 
on  the  principle  that  a  wave  incident  upon  a  boundary  is  partly  reflected 
and  partly  transmitted,  depending  on  the  boundary  conditions.  Previous 
applications  of  the  ultrasonic  reflection  approach  include  the  measure¬ 
ment  of  viscosity  in  a  liquid  using  obliquely  incident  shear  waves,!4)  and 
the  measurement  of  sound  velocity  in  solids  using  either  obliquely  inci¬ 
dent  longitudinal  wave  pulses (5)  or  normally  incident  longitudinal  contin¬ 
uous  waves.  (6) 

In  this  program  we  utilized  the  reflection  of  normally  incident  longitudinal 
wave  pulses  to  determine  fuel  density.  At  normal  incidence  the  reflection 
coefficient  for  plane  waves  depends  only  on  the  characteristic  acoustic 
impedances  of  the  adjacent  media  at  the  boundary.  These  media,  for 
example,  may  be  a  probe  and  the  fuel.  Since  the  fuel1  s  characteristic 
impedance  equals  pc,  if  the  probe'  s  impedance  Zj  is  known,  measure¬ 
ments  of  the  reflection  coefficient  at  the  boundary,  together  with  a  mea¬ 
surement  of  sound  speed  c  (or  time  interval  T,  where  T  is  proportional 
to  l/cl,  enable  one  to  determine  the  fuel  density  p. 

TASK  DESCRIPTIONS 


Having  generally  introduced  the  three  main  ultrasonic  approaches  used 
in  this  program  -  scattering,  transmission,  reflection  -  we  present  next 
a  brief  description  of  the  three  tasks  into  which  the  program  was  divided. 

Task  I  -  Theoretical  and  Analytical  Optimization  Study 

In  this  task  we  conducted  an  analysis  to  determine  the  optimum  parameter 
of  the  above  mass  flowmeter.  The  principal  factors  included  in  the  anaT 
y sis  were  temperature  range,  flow  range,  fractional  contaminant  content, 
pipe  geometry,  noise  level  as  a  function  of  frequency,  allowable  system 
response  time,  and  suitability  for  Army  aircraft  gas  turbine  installations. 

We  also  analyzed  six  piezoelectric  transducer  designs  as  candidates  for 
use  in  the  velocimeter.  They  are  listed  in  Table  II. 

The  normally  incident  longitudinal  wave  design  was  analyzed  for  the 
velocimeter,  to  convert  the  phase  shift,  which  is  proportional  to  v/ c^, 
to  flow  velocity  v.  This  design,  using  a  second  transducer  and  special 
probe,  was  also  analyzed  for  the  densitometer. 
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TABLE  II.  TRANSDUCER  APPROACHES 


Wave  Type 

Method  of  Introduction 

Longitudinal 

Phased  array 

Longitudinal 

Axially  incident  (from  end  wall) 

Longitudinal 

Obliquely  incident  (from  solid  wedge) 

Longitudinal 

Obliquely  incident  (from  recessed  liquid  wedge) 

Longitudinal 

Normally  incident  (from  side  wall) 

Shear 

Obliquely  incident  (from  solid  wedge) 

Task  11  -  Design  and  Fabrication 

This  task  involved  designing  and  fabricating  the  mass  flowmeter  system. 
The  mass  flowmeter  system  included  three  main  electronic  functions:  a 
velocimeter,  a  time  intervalometer ,  and  a  densitometer.  Additionally, 
several  transducers  and  measurement  cells  were  designed  and  built. 

Task  III  -  Tests 

In  this  task  we  tested  the  system  components  to  detect  and  correct  defi¬ 
ciencies,  to  determine  optimum  values  of  adjustable  operating  parameters, 
and  to  verify  performance.  We  also  conducted  room  temperature  tests  on 
'he  complete  system  at  flow  rates  of  50  to  5000  lb/ hr.  These  tests  in¬ 
cluded  the  use  of  water,  Avgas- 100,  JP-4  and  JP- 5  fuels,  and  a  fluid 
contaminated  as  specified  by  MIL-E-5007C.  The  purpose  of  these  tests 
was  to  determine  the  accuracy  of  the  mass  flowmeter  system  and  to  de¬ 
termine  the  effects  of  the  contaminant  on  the  mass  flowmeter  system. 
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ULTRASONIC  METHODS 


Following  the  format  of  the  previous  section  on  ultrasonic  approaches,  we 
describe  in  this  section  the  ultrasonic  methods  we  used  in  investigating 
scattering,  transmission  and  reflection.  These  methods  provided  the  test 
data  for  determining  flow  velocity,  density  and  mass  flow  rate.  A  dis¬ 
cussion  of  Reynolds  number  and  flow  profile  is  presented  first,  however, 
since  the  accuracy  of  any  ultrasonic  method  of  measuring  flow  velocity 
will  be  ultimately  limited  by  this  factor. 

REYNOLDS  NUMBER  AND  FLOW  PROFILE 


Reynolds  Number  is  the  dimensionless  ratio  Re  =  p  vD/n  =  vD/\<  .  For 
common  fuels,  the  temperature  dependence  of  p  and  v  are  shown  in 
Figures  1  and  2.  To  graphically  illustrate  the  change  in  Re  as  a  function 
of  fuel  type,  temperature,  and  flow  velocity,  the  Reynolds  number  nomo¬ 
gram  of  Figure  3  was  constructed.  The  extremes  of  Re  occur  when  v/v 
is  minimum  or  maximum,  as  shown  in  Table  III. 


TABLE  III.  RANGE  OF  REYNOLDS  NUMBERS  ENCOUNTERED 
IN  TWO  FUELS  FOR  M  FROM  50  TO  2000  LB/HR 
AND  TEMPERATURE  FROM  -65°F  TO  +160°F, 


CALCULATED  FOR  A  1/ 2-INCH  DIAMETER  PIPE 


_ 

Fuel  Type 

Minimum  Re 

Maximum  Re 

Re 

M 

Temp 

Re 

M 

T  emp 

(■) 

(lb/hr) 

(°F) 

(-) 

(lb/ hr) 

(°F) 

JP-4 

mm 

m 

•65 

2000 

160 

Avgas- 100 

mm 

mm 

-65 

2000 

160 

Experimentally,  the  flow  velocity  distribution  v(r)  is  observed  to  depend 
on  Re  as  shown  in  Figure  4.  The  turbulent  flow  data  may  be  represented 
by  an  empirical  equation  of  the  formal 

v(r)/v<0)  =  [(R-r)/R]l/n  (1) 

where  n  varies  slightly  with  Re.  n  =  6  at  Re  =  4000,  n  =  6,  6  at  Re  = 

23,  000,  and  n  =  7  at  1 10,  000.  The  i  atio  of  mean-to-maximum  velocity  is 

v  /  v ( 0 )  =  2n2/ (n  +  1)  (2n  +  1)  (2) 

d 
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KINEMATIC  VISCOSITY,  CENTISTOKES 


TEMPERATURE,  °C 


-60  -40  -20  0  20  40  60 


Figure  2.  Kinematic  Viscosity  (v  =  rj  / p  ) 
Versus  Temperature. 


STANCE 


For  n  =  6,  va/  v(0)  =  0,  791.  For  n  -  7,  va/  v(0)  =  0.  817.  For  laminar 
flow,  the  profile  is  parabolic: 

v(r // v (0)  =  1  -  (y/R  -  l)2-  1  -  (r/R)2  (3) 


and  v  / v(0)  =  0.  5. 
a 

If  one  interrogates  the  fluid  with  a  wave  that  effectively  samples  the  diam¬ 
eter,  but  not  the  whole  area,  one  measures  =  v  /  K,  where  K  is  less 
than  unity  and  depends  on  Re,  as  mentioned  before.  In  terms  of  the  pipe 
friction  coefficient  A  ,  K=  1/(1  +  0.  44  1 94  \rA).  For  A  between  0.  0  1  and 
0.  06,  K  ranges  between  0.  9023  and  0.  9577.  Various  oblique  paths  that 
measure  v^  are  shown  in  Figure  5,  and  their  projections  are  shown  in 
Figures  6  and  7.  Relative  to  a  clamp-on  flowmeter,  Figure  7d  shows  a 
v^  measurement. 


The  profile  correction  factor  K-  v^/  v^  may  be  expressed  in  terms  of  the 
velocity  distribution  v(r): 


K  - 


\  "R  -0 


L  /-R  -R 

~  /  2 it  rv(r )  dr  j  r 

R  Jq  2  •/) 


v(r ) dr 


if. 


R 


R  rR 


v(r)dr 


/ 
J  n 


v(r)  dr 


(4) 


For  parabolic  flow,  this  becomes 

/R  ,  ? 

r[  1  -  (r/R  ]  dr 

*  ±. _ - _ 

K  ”  R  /*R 


/ 

J  0 


=  3/4  =  0.  750 


[  1  -  (r/R)  ]  dr 


(5) 


while  for  turbulent  flow,  K  is  typically  about  0.  90  to  0.  96. 


Graphically,  deviations  due  to  the  dependence  of  K  on  Re  have  been  drawn 
by  Kritz  in  his  nomogram.  (&)  Kritz'  s  nomogram  is  derived  from  Prandtl's 
velocity- distribution  equation 


v(r) 


v(0) 


+  2.  5v* 


(6) 


where  v*  =  friction  velocity.  From  this, 

v  =  v  (1  +  0.  19  Re"°'  S 
d  a 


(7) 
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which  we  may  easily  relate  to  our 
1/(1  +  0.  19  Re"0,  1). 


K  -  v  / v  ,  if  we  write 

a'  d  Kritz 


Typical  errors,  according  to  Kritz,  are  listed  in  Table  IV,  together  with 
errors  we  calculated  for  laminar  flow.  The  discrepancies  for  laminar 
flow  may  be  explained  as  follows.  Prandtl1  s  equation  is  for  turbulent 
flow,  with  its  underlying  assumptions  most  valid  for  large  Re.(7)  Kritz 
erroneously  used  this  equation  for  Re  <  2000,  for  which  the  deviation  is 
in  fact  not  9.  5%  but  25%  ,  as  we  have  shown,  and  not  a  function  of  Re,  for 
Re  <  2000  (i.  e.  ,  laminar  flow).  Kritz1  s  nomogram  is  questionable  for 
transitional  flow  too  (2000  <  Re  <  40,000),  although  here  the  errors  would 
not  be  as  serious.  It  appears  valid  above  40,000,  except,  as  he  points  out, 
his  calculations  neglect  the  disturbance  to  the  effective  profile  due  to 
transducer  cavities  of  the  nonflush  type  (Figure  5a).  A  plot  of  the  error 
versus  Re,  for  Re  >  5000,  has  been  presented  recently  by  McShane,  (9) 
based  on  an  empirical  equation  of  the  form 


vd  =  v  (1.  1 19  -  0.  011  log  Re). 


(8) 


Table  IV  shows  that  errors  computed  from  (8)  or  (7)  agree  within  0.  5%  for 
104  <  Re  <  1CS. 


TABLE  IV.  CALCULATED  ERRORS  DUE  TO 


NONUNIFORM  FLOW  PROFILE 


Reynolds 

Number 

Kritz<8> 

Calculated  Error, 

McShane^ 

vd'va  (%) 

Present  Work 

1,  000 

9.  5 

25 

2,  000 

8.  9 

- 

25 

10, 000 

7.  6 

7.  5 

- 

50, 000 

6.  5 

6.  7 

- 

70,000 

6.  2 

6.  6 

- 

100,000 

6.  0 

6.  5 

- 

In  many  industrial  situations,  where  flow  is  turbulent,  temperature  and 
composition  of  the  process  fluid  are  sufficiently  constant  so  that  Reynolds 
number  can  be  estimated  from  va.  Since  K  is  relatively  insensitive  to  Re, 
the  estimated  Re  is  sufficient  to  determine  K  to  a  small  fraction  of  1%, 
which  means  v„  can  be  determined  to  that  accuracy  too. 

U  ' 


In  our  case,  however,  at  a  given  flow  velocity  v,  kinematic  viscosity  v 
and  hence  Re  can  change  by  an  order  to  magnitude  for  JP-4  over  the  full 


temperature  range  (Figures  1  and  2).  Therefore,  our  Re  compensator,  if 
used,  could  not  be  based  only  on  an  iteration  using  the  flow  velocity  v(r)  at 
a  point,  or  averaged  across  a  diameter,  if  errors  are  to  be  kept  below 
1%.  To  calculate  how  accurately  Re  must  be  determined  to  keep  errors  in 
K  less  than  a  prescribed  amount,  say,  0.  1%,  we  may  differentiate  the 
error  term  (1  -  v^/ va)  with  respect  to  Re.  This  leads  to  the  result  that, 
for  relatively  small  changes- in  Re,  the  percentage  change  in  K  is  approxi¬ 
mately  half  of  the  fractional  change  in  Re.  For  example,  if  ARe/ Re  =  0.  2, 
K  changes  by  about  0.  1%.  (For  an  order  of  magnitude  change  in  Re,  K 
changes  by  ~1%.  )  While  these  results  are  limited  to  turbulent  flow,  they 
suggest  that  determination  of  Re  to  _+  20%  could  permit  one  to  convert  to 
va  rather  accurately,  an  ultrasonic  measurement  of  v^. 

Let  us  examine  the  expression  Re  -  pv D/  rj .  Under  turbulent  flow,  v(0)  is 
equal  to  ~1.  16  v^  +  3%.  Assuming  p  and  D  are  known  to  better  than  1%, 
the  principal  uncertainty  in  Re  is  due  to  r) .  One  approach  to  obtaining  rj 
is  indicated  by  the  graph  in  Figure  8.  This  graph  is  constructed  using  the 
data  in  Figures  1  and  2,  and  shows  that  viscosity  may  be  computed  from 
measurements  of  density  and  temperature,  provided  the  fuel  behaves  as 
part  of  a  homologous  series.  (10) 

Another  approach  to  overcoming  the  uncertainty  in  Re  is  to  mix  the  profile 
to  achieve,  ideally,  K  =  1.  Assuming  one  uses  a  Kenics  mixer,  as  de¬ 
scribed  in  the  next  section,  one  can  calculate  the  distance  that  the  mixed 
or  flat,  uniform  profile  is  maintained.  For  low  flow  this  distance  is  about 
(.03)  (Re)  (diameter),  while  for  high  flow  thi3  distance  increases  to  ~50 
to  100  diameters. 

The  above  profile  considerations  may  be  summarized  as  follows: 

1.  Transducers  should  not  be  allowed  to  increase  the  uncertainty  in 
flow  profile  in  the  region  of  interaction  between  sound  wave  and 
flow. 

2.  For  the  range  of  conditions  expected  in  this  contract,  flow  profile 
could  require  compensation  (K)  factors  from  0.  75  to  0.  94. 

3.  For  laminar  flow  or  for  fully  developed  turbulent  flow,  but  not 
necessarily  for  transitional  flow,  the  correct  K  factor  could  be 
applied  to  the  extent  that  Re  were  known  and  to  the  extent  that  the 
profile  had  traveled  down  a  sufficient  length  of  smooth  pipe  to  be¬ 
come  established.  (As  a  test,  flow  could  be  measured,  in  principle, 
at  two  distances  down  the  pipe,  ) 

4  If  the  fuel  were  perfectly  mixed,  the  profile  would  be  uniform,  and 
K  would  iqual  unity.  This  condition  may  be  nearly  achieved  close 
to  the  exit  of  a  Kenics  static  mixer  at  high  flow,  and  to  a  lesser 
extent,  at  low  flow. 
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Figure  8.  Kinematic  Viscosity  Versus  Density-Temperature 
Function,  for  Six  Fuels. 
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5.  To  the  extent  that  scatterers  are  uniformly  distributed,  and 
uniformly  insonified,  the  average  Doppler  shift  can  be  processed 
to  yield  average  flow  velocity,  independent  of  the  flow  profile. 

6.  To  the  extent  that  a  plane  wave  can  interrogate  the  entire  cross 
section  of  flowing  fluid  without  becoming  curved,  its  transit  time 
can  be  processed  to  yield  average  flow  velocity,  independent  of 
the  flow  profile.  This  new  result  is  of  major  significance. 

SCATTERING 


The  Doppler  scattering  method,  automatically  weights  the  flow  profile  cor¬ 
rectly,  provided  insonification  is  uniform  and  provided  scatterers  are 
uniformly  distributed.  To  distribute  scatterers  as  uniformly  as  possible, 
and,  incidentally,  to  flatten  the  profile,  we  used  a  Kenics  static  mixer. 
This  mixer  is  a  helical  flow  divider  with  a  180  twist  per  section.  It  is 
virtually  clog-free.  When  k/[  =  2000  lb/hr,  seven  sections  typically 
introduce  a  pressure  drop  of  ?.bout  6  psi  in  a  l/ 2 -  inch- diameter  pipe,  and 
about  2  psi  in  a  3/ 4- inch- diameter  pipe. 

To  optimally  locate  the  interrogation  zone,  one  needs  to  know  how  far  the 
flat  profile  is  maintained  after  exiting  from  the  mixer.  This  distance 
depends  on  Reynolds  number,  which  in  this  work  ranges  from  about  100 
to  about  70,  000. 

For  high  Reynolds  numbers  (>50,000),  the  initially  flat  profile  transforms 
to  its  final  shape  after  ~50  to  100  diameters.  Therefore,  measurements 
close  enough  to  the  mixer  to  avoid  the  transformation,  i.  e.  ,  within  5  to  10 
diameters  after  the  mixer,  will  see  the  flat  profile.  This  location  contri¬ 
butes  to  highest  accuracy  at  the  highest  flow  rate. 

For  low  Reynolds  numbers  (<  1000),  the  initially  flat  profile  will  transform 
(relaminarize)  to  parabolic  in  a  distance  of  ~(0.  03)  (Re)  (diameter),  or  3 
diameters  at  Re  =  100,  (?)  This  means  that,  as  before,  to  avoid  the  trans¬ 
formed  profile,  one  should  interrogate  as  close  as  possible  to  the  mixer. 
But  even  within  ~1  or  2  diameters,  since  the  profile  is  already  curved  at 
minimum  flow  rate,  small  errors  accrue  to  the  extent  that  illumination  is 
nonuniform.  Note  that  the  scatterers  should  still  be  uniformly  distributed, 
even  though  flow  profile  is  parabolic.  Therefore,  uniform  illumination 
could  avoid  the  profile- dependent  errors, 

A  diagram  of  a  dual- frequency  Doppler  flow  velocimeter  is  shown  in 
Figure  9.  Since  the  basic  Doppler  shift  is  proportional  to  v/c,  an  input 
in  terms  of  fuel  type  and  temperature,  or  c  (where  c  depends^)  on  fuel 
type  and  temperature,  as  in  Figure  10)  is  required  if  v  is  one  of  the  de¬ 
sired  outputs.  Figure  1  1  indicates  the  theoretical  improvement  which 
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Figure  9.  Doppler  Flow  Velocimeter  Using  Two 
Transmitted  Frequencies. 


23 


SOUND  VELOCITY  ,  rVs 


TEMPERATURE,  °F 

-65  0  60  120  165 


Figure  10.  Sound  Velocity  Versus  Temperature 
for  JF-4  and  Avgas-100.  ^ 
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Figure  11.  Comparison  of  Uncertainties  Obtained  With 
Single-  and  Double -F requency  Doppler  Flow 
Velocity  Meters. 
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may  be  obtained  by  processing  Dopplers  front  two  different  transmitted 
frequencies  li  and  1',. 

Scattering  and  Frequency  Selection 

Several  design  calculations  were  done  in  order  to  more  precisely  deter¬ 
mine  tlie  operating  parameters  of  the  dual  frequency  flow  velocity  mea¬ 
suring  system,  First,  on  the  basis  of  the  contaminants  specified  in 
MIL-E-5007C,  the  scattering  levels  per  unit  volume  of  fuel  due  to  each 
of  the  various  contaminant  categories  were  determined.  Second,  the 
overall  scattering  cross  section  per  unit  volume  of  fuel  at  20  Mil/,  was 
computed  and  found  to  be  -0.  3  x  10'-'  cm2/  ce.  Assuming  a  pipe  diam¬ 
eter  of  1/2  in,  and  an  interrogation  length  of  1/  2  in. ,  I  he  volume  ol  fuel  from 
which  scattering  is  being  measured  is  about  3cm*.  This  will,  therefore, 
provide  an  overall  hackscatter  cross  section  of  about  10”  ^  ui/  Third, 
tin'  received  power  due  to  this  bulk  scattering  cross  section  was  com¬ 
puted,  and,  assuming  a  delivered  power  to  the  transducer  of  one  wait  per 
channel,  a  s igna  1  -  to- noise  ratio  of  40  dD  was  determined,  See  Appendix  I. 

A  number  of  trade-ofls  exist  in  the  selection  of  operating  frequency. 

The »i?  are  summarized  below. 

The  frequency  f  ot  the  wave  to  be  scattered  should  be  high  enough  to  avoid 
most  engine  noise  and  to  avoid  reverberation  in  the  pipe  wall  and  fluid. 

Vet  it  should  not  be  so  high  that  attenuation  is  excessive.  Extensive 
datadD  is  available  on  «  in  liquids  similar  to  aviation  fuels.  For  our 
temperature  range  the  measured  10  *7  a/ f2  is  less  than  100.  This  means 
at  f  10  MHz,  or mnx  <  { 1 0 “  1 7 )  (107)2  (100)  =  0.  1  Np/cm  ~  1  tIB/cm.  At 
20  MU/.,  a  m  a  K  =  4  dB/  cm.  Since  the  total  path  is  expected  to  be  to 
10  cm,  the  maximum  usable  frequency  should  be  selected  in  the  above 
rarge,  namely,  1?  to  20  MHz.  («  is  denoted  the  attenuation  coefficii  lit .  ) 

Again,  the  acoustic  attenuation  coefficient  in  the  medium  is  roughly  pro¬ 
portional  to  f2.  Therefore,  subject  to  geometrical  constraints,  too  high 
an  operating  frequency  range  will  give  rise  to  a  prohibitively  low  signal- 
to-noise  ratio.  On  the  other  hand,  the  overall  scattering  cross  section  of 
the  medium  increases  with  frequency,  and  the  statistics  of  the  measure¬ 
ment  improve  with  increasing  frequency.  Consideration  of  these  factors 
leads  to  the  choice  of  operating  frequencies  at  approximately  20  MHz. 

Response  T inn; 

For  performance  diagnostics  in  this  contract,  response  time  of  0.  5  sec  is 
desired.  For  control, <30  ms  is  the  objective. 
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Those  values  wort*  interpreted  as  the  time  it  takes  the  output  signal  Lo 
reach  1/e  of  its  final  value,  in  response  to  a  step  change  in  mass  flow 
rate  M.  It  appears  that  only  the  flow  velocity  (v)  part  of  the  system  need 
be  this  fast,  as  it  seems  unlikely  that  fuel  temperature  (density,  viscosity) 
could  change  so  quickly.  That  is,  temperature  functions  (density ,  vis¬ 
cosity)  can  be  monitored  with  longer  integration  times,  with  smoothed 
outputs  always  available  for  multiplication  with  v,  to  yield  Ivl.  This  sim¬ 
plifies  the  determination  of  M  without  sacrificing  the  essential  fast 
response  to  sudden  changes  in  flow  velocity. 

The  minimum  response  time  is  limited  in  one  respect  by  the  number  of 
scatterers  interrogated  in  the  integration  period.  This  depends  in  part 
on  the  gate  width,  which  may  be  ~20  to  30p  s  wide.  For  colinear  trans¬ 
ducers,  if  the  gate  is  too  wide,  the  profile  may  be  smeared  towards  the 
transducer  end  of  the  gate,  which  smears  the  Doppler  returns.  If  the  gate 
is  too  narrow,  longer  integration  times  are  required,  to  sample  statis¬ 
tically,  enough  scatterers.  A  gate  width  corresponding  to  How  along  an 
axial  path  equal  to  one  diameter,  for  1/2-in.  -diameter  pipe,  would  be 
about  25p  s.  For  transducers  oriented  at  a  right  angle  to  one  another,  the 
intersection  of  the  two  beam  patterns  defines  the  measuring  /.one,  even 
for  a  wide  gate.  This  permits  the  use  of  continuous  waves,  without  gating. 

T  ransducers 

The  simplest  piezoelectric  transducers  are  the  plane  wave  types.  We  also 
considered  annular  types  such  as  the  phased  array  and  the  related  dif¬ 
fraction  grating.  Annular  types  could  be  either  phased  array  or  obliquely 
incident  shear.  However,  it  turned  out  that  the  simplest  plane  wave 
offered  two  important  features.  First,  it  permitted  the  use  of  continuous 
waves.  Second,  by  locating  one  of  the  transducers  in  Figure  12  or  13 
adjacent  to  the  interaction  (scattering)  zone,  the  total  path  x  is  nearly 
halved.  Tins  means  attenuation  may  be  reduced  by  a  factor  of  nearly  o'* 
or  7.  3.  Conversely,  for  a  given  allowable  total  attenuation  vax,  fre¬ 
quency  f  can  be  increased  by  e  --  2.  7,  since  a  is  proportional  to  f  .  Again, 
higher  frequency  leads  to  scattering  contributions  by  smaller  particles, 
and  provides  better  statistical  averaging  over  more  particles. 

One  potential  drawback  of  the  right  angle  transducer  configuration  is  due 
to  the  changing  profile  between  the  wall  and  the  measuring  zone.  Therefore 
we  analyzed,  for  transient  conditions,  the  effect  of  different  average  flow 
velocities  along  different  paths  between  interaction  zone  and  receiver 
transducer.  It  turns  out  that  such  effects  are  negligibly  small.  This  re¬ 
sults  from  the  dual-frequency  Doppler  output  being  of  the  form 

I  =  f  (SL  S2  -  S2  S,)  dt  (4) 

Jo 
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Figure  12.  Cross  Section  Assembly  View  of  Laboratory 
Test  Fixture  (Top)  for  Doppler  Flow  Veloci- 
meter  Using  Flowing  Liquids,  and  also 
(Bottom)  Verification  of  Doppler  Effect  on  a 
Single  Scatterer  in  Water  Under  No-Flow 
Conditions. 
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3.  Demonstration  of  Doppler  Effect  Under 
No-Flow  Conditions,  Using  Gas  Bubbles 
as  Moving  Scattercr. 
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where  the  signals  S  are  of  the  form 


Sj  =  sin  [  Ujt  +  (v  +  c)] 


q  -  - 

i  a  (v  +  c) 


A 

d  (v  +  c)  _  ^1^  d  (v  +  c) 

dt  .  .2  dt 

(v  +  c) 


cos  [A] 


and  S2  and  are  similar,  with  subscript  2  instead  of  1.  (R  is  the  range.  ) 
For  transients  shorter  than  the  system  response  time, 


(v  +  c)j=. 


V T  ■ 

f  A  wR  .  f 

=  J  FT T)  sln  [ 


A  «  + 


A  qjR 
(v  +  c) 


1  d  (v  +  c 

J  (v  +  c) 


(v  +  c). 


which  is  negligibly  small,  since  the  flow  velocity  v  <<  c,  the  sound  speed. 
In  the  case  of  the  desired  signals,  however,  we  have  integrals  of  the  form 


r 

=  / 


(v  +  c) 


sin  A  wt  + 


A 

,  u>R  j 

+  c)  J 


which  are  relatively  large  since  R(T)-R(0)  is  substantial.  The  system  out¬ 
put  would  thus  depend  mainly  on  flow  in  the  selected  scattering  zone,  and 
not  on  flow  velocity  fluctuations  outside  this  zone. 


Scattering  Observations 


The  method  and  equipment  used  to  observe  scattering  from  moving  par¬ 
ticles  are  shown  in  Figures  12  and  13.  The  oscillator  provides  a  stable 
source  of  continuous  waves  essentially  independent  of  transducer  tempera¬ 
ture.  Laboratory  tests  were  conducted  on  moving  scatterers  in  a  station¬ 
ary  liquid,  and  also  on  scatterers  carried  at  the  moving  iiquid'  s  flow 
velocity.  These  tests  used  uncontaminated  Avgas- 100  and  also  water  as 
the  fluid  media.  In  these  tests,  using  20  MHz  continuous  waves,  scattered 
waves  were  not  detected  except  when  one  or  more  bubbles  were  present 
naturally  or  present  by  controlled  introduction.  Calculations  predicted 
that  fuel  contaminated  above  1%  of  the  maximum  levels  specified  in  MIL- 
E-5007C,  paragraph  3.  4,  1.  3,  Table  1,  would  yield  detectable  scatter. 
However,  the  need  for  high  accuracy  and,  ultimately,  fast  response, 
means  that  scatter  must  be  not  merely  detectable,  but  ample.  Instead  of 
pursuing  methods  of  increasing  the  scatter,  or  increasing  sensitivity  to 
scatter,  we  reexamined  alternative  transmission  approaches, 
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TRANSMISSION  METHODS 


Various  ultrasonic  transmission  methods  were  analyzed  for  measuring  v 
and  c.  In  many  cases  one  would  prefer  to  either  eliminate  the  need  for 
a  c  measurement,  or  choose  a  method  that  yields  both  v  and  c.  In  this 
program,  however,  our  approach  required  both  variables  to  be  deter¬ 
mined  independently. 

Flow  Velocimeter 


Translation  of  the  40:1  M  range  of  50  to  2000  lb/ hr,  for  fuels  having 
densities  ranging  from  about  0,  65  to  nearly  0.  90  g/ cn/,  leads  to  the 
requirement  to  measure  flow  velocity  over  a  ~60:1  range,  from  ~0.  05 
to  m/  s  (~2  to  ~120  in./  s)  ir.  a  1.  27-cm  (1/ 2-inch)  diameter  pipe.  A 
nomogram  relating  M,  v,  p  ,  pipe  diameter,  and  volumetric  flow  rate  ai 
2000  lb/  hr  is  shown  as  Figure  14. 

Pulse  Time  of  Flight 

Let  us  first  estimate  the  axial  path  length  L  required  to  meet  accuracy 
requirements  of  5%  at  minimum  flow  and  0.  5%  at  maximum  flow.  The 
upstream  and  downstream  travel  times  are  L/  (c  4  v)  and  L/  (c  -  v). 
Their  difference  is  At  ~  2Lv/  c^,  from  which  L  =  At/2v.  Assuming 
a  time  interval  resolution  of  +  1  ns,  and  c  =  1500  m/ s,  we  compute 
L  =  45  cm  (~18  in.  )  when  v  =  ,  05  m/  s,  and  L  =  7.  5  cm  ( ~3  in.  )  when 
v  =  3  m/  s.  These  L'  s  are  minima,  since  they  cause  us  to  use  up  our 
entire  M  error  allotment  in  determining  v,  without  allowing  any  mar¬ 
gins  of  error  in  p  or  c. 

If  we  allocate  half  our  M  error  to  v,  then  we  must  double  L  to  90  cm 
(~36  in.  )  to  meet  the  accuracy  requirement  at  minimum  flow.  If  the 
unwieldy  values  of  L  between  45  and  90  cm  are  not  yet  sufficient  reason 
to  force  abandoning  this  method,  two  further  considerations  certainly 
clinch  this  decision.  First,  to  achieve  _+  1  ns  resolution,  the  pulse 
width  should  be  100  ns  or  shorter,  corresponding  to  a  frequency  of 
M0  MHz.  Attenuation  of  10  MHz  pulses  is  estimated  to  be  M  dB/ cm, 
which  value  leads  to  signal  loss  of  45  to  90  dB  in  the  fuel. 

Second,  in  propagating  over  such  long  paths,  well  into  the  far  field 
(Fraunhoffer  zone)  the  initially  plane  wavefront  will  be  increasingly 
distorted  at  increasing  velocities  of  laminar  flow  (Re — *2000)  and 
somewhat  distorted  even  in  the  flatter  profile  associated  with  turbu¬ 
lent  flow  (Re  —•70,  000).  These  considerations  show  that  while  one 
might  achieve  acceptable  accuracy  over  the  upper  half  of  the  flow 
range  using  an  L  <  20  cm  (~8  in.  ),  the  pulse  time  of  flight  method 
is  unsuitable  for  the  lower  flow  velocities.  Appendix  II  supports  thi a 
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conclusion  with  an  analysis  of  far  field  errors  under  laminar  flow 
conditions. 

Since  wavefront  curvature  under  laminar  flow  forces  us  to  abandon 
long  paths,  and  the  accuracy  requirement  at  low  flow  forces  us  to 
abandon  the  pulse  time  of  flight  method  for  short  paths,  we  are  appar¬ 
ently  left  with  only  two  alternative  short  path  ultrasonic  transmission 
methods:  phase  measurement  at  fixed  frequency,  or  frequency  mea¬ 
surement  at  fixed  phase.  In  the  present  program,  we  investigated  the 
former  alternative. 

Continuous  Wave  Phase  Measurement 


Based  on  the  values  of  v,  c,  a  ,  taking  10  cm  (~4  in.  )  as  the  maximum 
cell  dimension,  and  setting  the  maximum  phase  shift  difference  equal 
to  180  ,  we  selected  5  MHz  as  the  operating  frequency.  Measurement 
of  this  phase  difference  to  +.  0.  3  corresponds  to  an  error  in  vmax/  c^ 
of  about  0.  2%.  At  the  low  flow  limit,  _+  0.  3°  corresponds  to  an  error 
of  about  _+  10%.  Phase  meters  are  now  available  with  digital  and  analog 
phase  accuracy  of  0.  3°  at  0.  0  1  to  3V  single  frequencies  up  to  1  1  MHz 
(e.  g.  ,  Dranetz  Engineering  Laboratories'  model  305C  phase  meter 
with  305-PA-3005  plug-in).  It  therefore  appears  that  errors  in  flow 
velocity  of  0.  5%  or  less  at  vmax,  and  10%  at  vmjn,  are  to  be  associated 
with  presently  available  electronic  instruments  that  measure  phase  shift. 
By  doubling  the  operating  frequency  to  10  MHz,  to  provide  360°  phase 
difference  at  vmax  and  b°  at  vmjn,  it  appears  that  these  percentage 
errors  could  be  halved,  to  0.  25%  and  5%,  respectively, 

The  measurement  to  be  described  next  involves  the  processing  of  essen¬ 
tially  plane  5  MHz  continuous  waves  (cw)  propagated  in  the  transducer1  s 
near  field  through  100%  of  the  flow  cross  section.  Flow  is  measured  in 
a  square  channel,  with  a  rectangularly-collimated  beam  inclined  at  45° 
to  the  flow.  This  approach,  while  it  may  still  need  some  refinement 
relative  to  inlet  geometry  and  channel  interruptions,  provides  a  new 
instrumentation  method  for  responding  essentially  linearly  and  accu¬ 
rately  to  flow  velocity  in  transmission  measurements  over  laminar, 
transitional  and  turbulent  flow  conditions,  without  requiring  knowledge 
of  or  compensation  for  the  Reynolds  number  or  flow  profile. 

Similar  to  previous  ultrasonic  transmission  flowmeters,  the  system 
to  be  described  now  measures  the  phase  difference  between  waves  trans¬ 
mitted  upstream  and  downstream.  The  two  cw  waves  traverse  virtually 
the  same  path  at  virtually  the  same  time.  The  two  waves  have  the  same 
carrier  frequency,  yet  can  be  separated,  or  identified,  by  the  manner 
in  which  they  are  each  modulated  and  then  demodulated. 
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This  advanced  flow  velocimeter  is  shown  in  Figure  15.  In  this  system, 
a  5  MHz  crystal  controlled  oscillator  signal  is  divided  into  two  isolated 
channels.  The  waves  in  the  upper  and  lower  channels  are  phase  coded 
by  means  of  a  dual  coherent  pseudo- random-noise  code  generator, 
producing  a  sequence  of  l's  and  0's  such  that  the  resulting  code  dis¬ 
plays  all  the  characteristics  of  perfect  randomness.  The  length  of  the 
code,  as  well  as  the  clock  rate  (which  is  derived  by  dividing  down  from 
a  6.  2  MHz  carrier)  are  arbitrarily  variable.  Each  "1"  produces  a 
+  90°  phase  shift  of  the  carrier,  while  each  "0"  produces  a  -90  phase 
shift.  The  resulting  spectrum  strongly  resembles  a  gaussian  noise 
spectrum  centered  at  the  carrier  frequency  and  has  a  half  width  equal 
to  the  reciprocal  of  the  clock  rate.  The  codes  in  the  upper  and  lower 
channels  are  different  codes.  These  pseudo- random  sequences  have 
the  property  of  being  orthogonal,  which  means  that  a  signal  coded  with 
code  Cj  but  detected  with  a  code  C£  will  be  rejected  by  a  factor  1/N  in 
voltage  (where  N  =  number  of  bits  in  the  code)  as  compared  with  a  sig¬ 
nal  coded  with  Cj  and  detected  by  Cj.  Thus,  a  code  having  10^  bits 
will  have  a  correlation  function  of  40  dB;  while  10^  bits  will  provide  a 
correlation  function  of  80  dB!  '"  Another  property  of  these  (shift  reg¬ 
ister)  codes  is  that  a  code  which  is  shifted  by  one  or  more  bits  is 
orthogonal  to  the  unshifted  code.  Thus,  it  is  possible  to  provide  range 
gating  of  a  high  order  of  isolation,  as  well  as  the  isolation  of  two  sep¬ 
arately  coded  signals.  The  use  of  phase  coding  permits  simultaneous 
and  continuous  operation  of  one  or  more  systems  operating  on  exactly 
the  same  basic  carrier. 

Going  back  to  Figure  15,  we  see  that  the  two  5  MHz  waves  are  indepen¬ 
dently  phase  code  modulated  (C^  and  Cg)  by  the  double  balanced  mixers 
(DBM1  s)  at  the  left.  These  are  then  fed,  via  circulators  (which  are 
actually  hybrid  power  dividers)  to  the  two  transducers.  The  received 
waves,  which  of  course  resemble  wide  band  noise,  are  fed  (via  the 
circulators)  to  the  top  and  bottom  double  balanced  mixers  on  the  right. 
The  demodulating  codes,  and  Cj^,  which  are  delayed  just  the 

proper  amount,  are  also  fed  to  these  DBM1  s.  The  outputs  of  the  DBM's 
are  two  cw  signals,  as  well  as  wide  band  noise.  The  "noise"  results 
from  noncorrelation  of  unw  .ted  signals  from  either  the  wrong  channel 
or  from  stray  reflections  and  r  ver berations  within  the  pipe.  Only  that 
wave  which  is  due  to  a  reverbere  io n-  fte  str aight-through-transmission 
results  in  a  cw  wave  at  the  dome  ul  tor  DBM' s.  We  ne::t  go  through  a 
pair  of  tuned  5  MHz  amplifiers  which  stri”  off  most  of  the  uncorrelated 
information,  and  finally  feed  into  a  pair  of  phase-locked  loops.  These 
devices  are  essentially  se  f-tracking  filters  and  employ  a  voltage  con¬ 
trolled  oscillator  (VCO)  ana  a  phase  detector  in  a  feedback  configuration 
such  that  the  VCO  is  forced  to  track  the  cw  input  signal.  The  bandwidth 

*An  electromagnetic  system  recently  completed,  using  coding  tech¬ 
niques  identical  to  this,  exhibits  80  dB  correiation/decorrelaticn  ratio. 
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Figure  1  5. 

Block  Diagram  of  Ultrasonic  Mass  Flowmeter. 

of  the  feedback  loop  is  adjustable  to  be  arbitrarily  narrow,  e.  g.  ,  100 
Hz.  The  outputs  of  the  phase-locked  loops  are  square  waves  whose 
phases  are  equal  to  the  phases  of  the  cw  +_  up  to  100  Hz  components  of 
the  input  signals,  and  whose  amplitudes  are  independent  of  the  input 
amplitudes.  These  square  waves  are  next  multiplied  in  the  right-hand 
DBM  to  produce  a  voltage  which  is  precisely  proportional  to  the  phase 
difference.  Instead  of  phase-locked  loops,  one  may  use  a  standard 
phase  meter.  Another  simplification  involves  using  just  one  code  for 
both  directions. 

One  might  ask  about  the  residual  uncorrelated  signal  within  the  100  Hz 
bandwidth  of  the  phase-locked  loops.  The  answer  is  that  there  isn't 
any!  The  reason  is  that  the  overall  time  duration  of  the  sequence  is 
made  to  be  short  enough  so  that  the  frequency  components  associated 
with  uncor related  signals  fall  well  outside  the  bandpass  of  the  feed¬ 
back  loop  in  the  phase-locked  loops. 

.  2 

To  obtain  M,  the  v/c  output  of  the  flow  velocimeter  is  combined  with 
T  and  p  c  as  shown  in  the  lower  portion  of  Figure  15.  Details  on  the 
T  and  p  c  measurements  are  given  below. 


Time  Intervalometer 


Measurement  of  the  time  interval  T  across  a  fixed  or  known  path  in  a 
direction  orthogonal  to  the  flow  is  required  in  order  to  obtain  M  from 
measurements  of  v/c2  and  pc.  That  is,  M  is  proportional  to  (l/T) 

(v/  c2)  (pc),  or  M  =  kVp/ cT  where  k'depends  on  the  pipe  area  and  units 
of  measure. 


The  method  used  to  measure  T  is  illustrated  in  Figure  16.  It  utilizes 
a  damped  5  MHz  transducer  mounted  outside  the  flow  channel.  The  trans¬ 
ducer  is  energized,  and  selected  echoes  of  like  polarity  are  measured, 
with  a  modified  Panametrics  5220  gage  operated  as  a  time  intervalometer. 
This  gage  automatically  averages  the  transit  time  between  ten  pairs  of 
selected  echoes.  The  reason  the  echoes  need  to  be  selected, is  that  follow¬ 
ing  the  initial  energizing  pulse,  the  wall  supports  many  reverberations  of 
the  pulse.  The  frequency  is  chosen  high  enough  so  these  are  substantially 
attenuated  by  the  time  echoes  arrive  from  the  far  wall  of  the  channel. 

Since  the  5220  is  designed  to  read  between  echoes  of  only  positive  polarity, 
the  measurement  is  made  between  two  successive  reverberations  in  the 
fluid,  rather  than  between  the  first  pair  of  wall/ fuel,  fuel/ wall  echoes. 

If  the  transducer  were  in  direct  contact  with  the  fuel,  T  could  be  measured 
from  the  initial  pulse  to  the  first  echo  from  the  fuel/ wall  interface. 
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REFLECTION  METHODS 


The  present  ultrasonic  reflection  coefficient  densitometer,  where  R  = 

(p  2C2  ~  ^  l )/ (P  2C2  *  ^0’  may  ke  comPared  with  two  other  types  of  ultra¬ 
sonic  densitometers.  In  one  type,  first  investigated  in  connection  with 
mass  flowmeters  by  Kritz(^)  and  later  by  Dalke  and  W'elkowitz, U2)  the 
output  voltage  in  a  resonant  circuit  is  a  measure  of  the  liquid1  s  charac¬ 
teristic  impedance  p  2C2*  This  p  ^2  *s  couP^ec^  into  the  resonant  circuit 
by  its  loading  effect  on  a  piezoelectric  crystal  which  radiates  into  the 
liquid.  For  operation  at  high  pressure,  a  half-wave  protection  plate 
separates  the  crystal  from  the  liquid. 

In  a  second  type,  the  distributed- spring-mass  vibration  densitometer, 
now  commercially  available  from  ITT- Barton  (Series  650),  the  resonant 
frequency  of  a  metal  tube  is  measured,  the  square  of  this  frequency  being 
inversely  proportional  to  the  fluid  density  p  3-  As  presently  constructed, 
this  tube  of  monel  and  316  SS  is  normal!/  centered  within  a  pipe  tee 
section  of  3.  8  cm  (1.  5  in.  )  diameter.  Since  the  tube  materials  have  moduli 
of  elasticity  which  are  temperature- dependent  (readout  could  vary  up  to 
_+  1%  over  range  -65°to  +  160  F),  some  correction  is  required  for  operation 
over  the  full  temperature  range.  (One  approach  to  minimizing  the  tem¬ 
perature  dependence  of  a  resonating  probe  is  to  build  the  probe1  s  reso¬ 
nant  shell  out  of  a  quartz  crysta1  whose  longitudinal  axis  is  inclined  at 
about  5  to  the  Z-axis,  or  optical  axis,  of  the  crystal.  (13)  )  Provided  the 
ITT  densitometer1  s  size,  intrusiveness  and  temperature  -  dependence  are 
allowable,  however,  its  precision  of  +0.001  specific  gravity  units  for  liquids 
is  attractive.  Effects  of  vibration  and  contaminants  should  be  negligible. 

In  contrast  to  these  two  types,  the  reflection  coefficient  densitometer 
probe  can  respond  to  a  very  small  volume  of  liquid  (e.  g.  ,  using  focussed 
waves,  it  could  theoretically  measure  a  droplet);  it  can  operate  over  wide 
pressure  ranges  without  compensation,  and  when  the  probe  is  fabricated 
of  certain  materials  such  as  AT  quartz,  it  can  operate  over  wide  tem¬ 
perature  ranges  without  compensation.  It  can  be  recessed  or  fiush- 
mounted.  Sensitivity  is  limited  by  echo  amplitude  resolution;  using  5  MHz 
video  pulses,  density  should  be  resolvable  to  about  0,  2%, 

Referring  to  Figure  17,  we  show  the  reflection  coefficient  densitometer 
schematically  in  two  configurations.  (Figure  16  shows  the  actual  cross 
section  as  used  in  the  mass  flowmeter.  )  The  stepped- diameter  probe  yields 
echoes  A  and  B  which  ring  for  one  or  more  cycles,  depending  on  acoustical 
and  electrical  terminations.  A  Panametrics  pulser/ receiver,  model 
5050  PR,  was  used  to  obtain  these  echoes.  In  Figure  16  the  time  interval 
across  the  liquid  is  measured  separately.  This  is  accomplished  by 
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Figure  17.  Schematic  of  Ultrasonic  Densitometer  and 
Sound  Velocimeter  Measuring  Cells. 
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using  the  lower  iriiHmlm  rr  .mil  it  modified  Danamel  i* if  *  thickness  gage, 
mode  I  to  measure  between  two  tie  let  l  i'll  eehneS  ol  like  polarity,  to 

j_  10  nu,  ( I •*) 

The  reflection  part  ut’  the  density  measurement  responds  to  p  >c*/Zj, 

This  i  product  in  collet, I  the  vhacuelert  t*lic  impedance  of  the  (fuel) 
medium.  Ity  measuring  the  amplitude  of  the  wave  reflected  >»t  the  pipe/ 
fuel  interfiice,  p^e,  is  determined  in  terms  of  Zj,  To  the  extent  that 
sensed  neitr  the  pipe  wall  in  r epresentative  of  t h **  average  in  the  fluid, 
the  method  should  provide  a  valid  measure  oi  fuel  density  (Figures  Ut,  1?), 

The  t rnnsmissiun  part  of  the  density  nu. inurement  responds  to  e^  alone, 
by  arranging  to  measure  the  time  T  tt  takes  for  a  sound  puUe  to  traverse 
the  test  cell,  Fur  a  measurement  aeross  the  inside  diameter  1),  the  sound 
speed  is  c  l  D/T,  For  simplieity.  let  us  assume  A  s  11  when  p  0, 

To  cunt  put  e  fuel  density  we  simply  take  tin*  ratio 

P  ,  1  /.A,)  (I  +  H)/tl  -  K)  (/./c,)l(A  -  M)/(A  +  H)  ]  (Id) 

ts  I  fc*  l  U 

where  II  is  tin'  sotuul  pressure  reflection  coefficient  at  the  prohe/fnel 
interface,  Calculations  show  that,  for  interface  reflections  of  ..uoul  IV, 

A  and  H  echo  amplitttde  measurements  accurate  tut),  l  ,uV  will  determine 
the  reflection  coefficient  It  to  0,  t"‘,i,  ami  p  a',  to  0,  ri"'„,  These  calculations 
assume  an  all*  stainless  -  steel  pipe  or  test  call  and  plane  waves.  One  can 
improve  on  the  p^c,  determination  by  us  inn  a  probe  of  characteristic  int- 
pedant  e  closer  to  that  of  the  fuel.  For  example,  using  a  fused  silica,  AT 
quart's,  ur  T-dO  glass  probe,  p  ^c^  and  hence  p  improves  to  M),  V*', ,  For  a 
plastic  probe  one  could  resolve  p  ^  to  0,  1%,  theoretically,  using  this  re- 
Hectometer  method. 

A  more  detailed  analysis  of  this  method  is  given  below, 


Ultrasonic.  Densitometer  fsenailiviti 


The  sound  pressure  reflection  coefficient  If.  at  the  buffer/ fuel  interface  is 
11  -  (r-l)/(r  +  l)  where  r  impedance  ratio  Z^/Zj,  7^  is  the  fuel  charac¬ 
teristic  impedance  p^c^,  and  7.  j  is  the  buffer  characteristic  impedance 
p  |t  |.  What  value  of  7.  |  maximizes  the  sensitivity  of  this  measurement 
to  fuel  density  ?  That  is,  what  value  oi  7.  j  maximizes  .sensitivity  with 
respect  to  7.  or  r? 

i* 


Since  we  propose  to  measuir  R,one  may  first  suppose  that  maximum 
sensitivity  occurs  when  dll/  dr  iu  maximum  (see  Figure  1H).  This  maximum 
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Figure  18,  Reflection  Coefficient  Functions. 


is  readily  found  to  occur  at  r  =  0,  for  which  he  slope  of  R  versus  r  is 
dR/  dr  -■  2/ (r  +  1)2  =  2.  However,  near  r  =  0,  R  =  -1  While  AR  is 
indeed  largest  here,  for  a  given  Ar,  we  would  be  sensing  a  relatively 
small  change  AR  in  a  large  number  R. 

Thus  we  are  led  to  consider  (l/R){dR/dr)  as  the  densitometer  sensitivity 
function  to  be  maximized  with  respect  to  r.  This  has  maxima  when  r=  1, 
that  is,  when  Z ^  s  Z ^  This  function  has  positive  and  negative  branches, 
as  is  readily  seen  by  evaluating  (l/R)(dR/ dr)  -  2/  (r^  -  l)  for  r  over  a 
range  straddling  r  =  1  (see  Figure  19).  In  this  work  we  are  concerned 
with  the  negative  branch,  since  fuel  impedances  are  less  than  solid 
buffer  impedances  Zjj  i.  e.  ,  r  =  Z^/Zj  <  l. 

The  above  analysis  leads  us  to  conclude  that  as  far  as  ultrasonic  reflection 
coefficient  measurements  of  fuel  density  are  concerned,  the  buffer  rod 
should  preferably  be  chosen  to  have  an  impedance  as  close  as  practical 
to  the  fuel  Z^. 

In  any  given  installation,  other  requirements  may  dictate  how  close  the 
design  can  approach  this  "acoustic"  objective. 

For  comparison  purposes,  we  list  the  characteristic  impedances  of 
Avgas- 100,  JP-4,  water  and  several  solids,  at  room  temperature,  in 
Table  V. 


TABLE  V.  CHARACTERISTIC  IMPEDANCES 

Medium 

Characteristic  Impedance, 
g/  cm2-p  s 

Avgas- 100 

0.  087 

JP-4 

0.  095 

Water 

0.  150 

Polystyrene 

0.  25 

Acrylic  (Lucite,  Plexiglas) 

0.  32 

Teflon 

0.  59 

Magnesium 

1.  01 

Pyrex 

1.  24 

Fused  Silica 

1.  30 

Aluminum 

1.  75 

T  itanium 

2.  77 

Stainless  Steei  302 

4.  55 

Tungsten 

9.  98 
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We  describe  next  the  main  considerations  underlying  the  optimization  of 
the  reflection  method. 

To  select  the  proper  material  for  the  probe,  consider  the  following:  A 
plastic  buffer  (Zj  =  0.  3  g/cm^-ps)  provides  relatively  good  sensitivity, 
but  unfortunately  plastics  are  too  temperature- sensitive.  To  temperature- 
compensate  to  _+  0.  1%  for  a  plastic  such  as  delrin,  one  would  need  to 
monitor  the  probe  temperature  to  about  1°F. 

A  better  approach  is  to  use  a  more  brittle  material  such  as  glassy  carbon 
(Zj  =  0.  7)  or  fused  silica  (Zj  =  1.  3)  whose  Zj1  s  change  very  little  fiom 
-65°  to  +160  F.  Z1  in  these  buffer  materials  is  1  to  2  orders  of  magnitude 
less  sensitive  to  temperature.  The  trade-off  is  a  factor  of  2  or  4  less  R 
sensitivity  to  Z^.  Fused  silica  is  known  to  exhibit  a  small  but  ultimately 
undesirable  temperature- dependence,  since  its  sound  velocity  Vj^  increases 
with  temperature  at  the  rate  of  ~8  1  ppm  per  °C.  Over  a  temperature 
range  of  ~120  C  (-65°  to  +160  F),  this  leads  toAV^/Vj^-  1%.  To  avoid 
the  need  for  electronic  temperature  compensation,  we  fabricated  identically- 
shaped  probes  of  AT  quartz  and  T-40  glass.  A  comparison  of  these  mate¬ 
rials  is  given  in  Table  VI. 

The  AT  quartz  has  the  lowest  AVj^/ AT,  but  must  be  carefully  oriented. 
Therefore  the  T-40  glass,  which  is  isotropic,  may  provide,  ultimately, 
a  more  practical,  and  certainly  lower  cost,  solution. 

Zj  and  C2  can  be  measured  with  errors  in  the  range  of  0.  1  to  0.  01%.  There¬ 
fore  the  practical  limit  in  this  method  is  imposed  by  one'  &  ability  to  peak- 
detect  and  measure  A  and  B,  or  their  ratio,  A/B.  If  A  and  B  are  video 
(nonoscillatory)  shaped,  peak  detector  accuracy  of  better  than  0.  5%  is 
achievable  at  present  only  if  A  and  B  are  at  least  1  pc  s  wide.  Such  wide 
pulses  can  be  generated,  but  they  do  not  propagate  as  plane  waves,  and 
therefore  may  complicate  the  probe  design.  However,  if  A  and  B  pulse 
shapes  remain  constant,  a  circuit  that  measures  a  fixed  proportion  of  the 
peaks  will  be  adequate  (Figure  20).  Also,  if  A  and  B  were  rf  bursts,  say, 

10  cycles  of  a  10  MHz  carrier,  the  detected  1  ps  envelopes  can  be  mea¬ 
sured  by  1  MHz  peak  detectors,  or  by  a  standard  instrument  such  as  a 
Matec  2470  Attenuation  Comparator,  which  resolves  0.  01  dB  v/ith  a  2,  5s 
response  time. 

The  circuit  of  Figure  20  operates  as  follows.  On  receipt  of  a  plus  or  minus 
trigger,  a  bipolar  switch- selectable  comparator  starts  a  timing  sequence. 
This  sequence  is  adjusted  so  that  two  2  ps  gates  straddle  the  A  and  B 
echoes  from  the  densitometer  probe.  The  peak  detectors,  built  using 
Fairchild1  s  pA7l0  and  Dynamic  Measurement'  s  FST-160B  generate  dc 
voltages  proportional  to  A  and  B  peak  magnitudes, 
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TABLE  VI.  COMPARISON  OF  DENSITOMETER  PROBE 
MATERIALS:  FUSED  QUARTZ,  AT  QUARTZ, 

T -40  GLASS 

Material 

Velocity 
Temp.  Coef. 
(ppm/°C) 

Expansion 

Coefficient 

(ppm/oQ 

Density 
(g/  cm3) 

Longitudinal 
Velocity 
(cm/p  s) 

Char.  Acoust. 
Impedance 
(g/cm2-ps) 

Fused  Quartz* 

81.  4 

0.  55 

2.  20 

0.  590 

1.  30 

AT  Quartz** 

-2.  0 

a  3=  7.  8 
a  j  =  a  ->  =  14.  3 

2.  65 

0.  688 

1.  85 

T-40  Glass*** 

-2.  8 

~8 

3.  39 

0.  435 

1.  47 

*Data  from  Amersil,  Inc.  ,  685  Ramsey  Avenue,  Hillside  N.  J. 

**Data  from  W.  P.  Mason,  Piezoelectric  Crystals  and  Their  Application 
to  Ultrasonics,  102-105,  Van  Nostrand  (1950),  longitudinal  mode. 
vv>.<Data  from  Bausch  &  Lomb,  Inc.  ,  Special  Products  Division,  635  St. 
Paul  Street,  Rochester,  N.  Y.  14602;  private  communications , 

Dr.  Hensler.  Average  temperature  coefficient  of  -2.  8  ppm/°C  is 
over  the  range  -20° to +37  C,  with  zero  point  of  37  C.  Zero  point 
temperatures  vary  from  -10°C  to  70°C.  See  also:  J.  T.  Krause, 

IEEE  Trans.  Sonics  and  Ultrasonics  SU-  19  (1),  34-40  (Jan.  1972); 

J.  T,  Krause,  private  communication. 


The  sum  and  difference  amplifiers,  National  301  A'  s,  together  with  a  "k" 
potentiometer,  yield  the  terms  kA  _+  B.  The  quotient  of  these  terms,  when 
multiplied  by  the  time  intervalometer  output  T,  yields  a  result  proportional 
to  density  p  An  operational  amplifier  scales  this  result,  to  end  up  with 
the  final  output  value  for  p  2 . 

The  circuit  is  aligned  using  the  following  laboratory  procedure, 

1)  Select  proper  sync  polarity  with  switch  located  at  back  of  board. 

2)  Set  Trigger  Level  potentiometer  to  reject  any  noise  in  sync  signal 
if  noise  is  present. 

3)  Set  Delay  potentiometer  so  that  gate  pulses  at  TP6andTF7  overlap 
display  of  desired  portions  of  waveform. 

4)  Turn  off  source, 

5)  Set  kA  to  0  (full  ccw). 
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6)  Set  kA  +  B  offset  potentiometer  for  0  mV  at  TP9  (Pin  8). 

7)  Set  kA  -  B  offset  potentiometer  for  0  mV  at  TP10  (Pin  9). 

8)  Turn  source  on  and  adjust  kA  potentiometer  for  null  at  TP10. 


Densitometer  Noise  Analysis 


If  noise  were  absent,  the  density  would  be  determined  as 


But  suppose  noise  of  magnitude  e  is  added  to  A  and  subtracted  from  B. 
Then 


P2U)  =  — 


(A  +  g  )  -  (B  -  e  ) 
(A  +  c)  +  (B  -  e) 


A  -  B  -  2t 
A  +  B 


(17) 


Let  tlL  .  P2(€)~P2  -2e 

p  2  ”  p  2  A-B 

Now  if  p^  <  0.  2%,  and  if  A  =  IV,  B  »  0.  9V,  we  require  that 

t  <  100  p  V. 

This  analysis  shows  that  error  voltages  up  to  0.  1%  of  (A-B)  are 
tolerable. 


(18) 

(19) 


OTHER  ULTRASONIC  METHODS 

The  foregoing  discussion  has  covered  those  scattering,  transmission  and 
reflection  methods  which  comprised  the  major  effort  in  this  program.  It 
is  to  be  understood,  however,  that  several  variations  may  be  of  interest 
in  the  future,  depending  on  changes  in  technical  requirements  regarding 
the  types  of  fuels,  flow  range,  temperature  range,  accuracy,  response 
time,  geometry,  etc.  Such  changes  might  occur,  for  example,  as  a  re¬ 
sult  of  emphasizing  equipment  simplicity  rather  than  emphasizing 
versatility. 

To  illustrate  the  nature  of  these  variations,  we  briefly  indicate  several 
other  ultrasonic  methods  of  measuring  v  or  v/c^,  and  p  . 
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Determination  of  Flow  Velocity  v 


A  direct  determination  of  v,  without  requiring  a  c  measurement,  may  be 
achieved  by  controlling  two  variable  frequencies  f^,  such  that  the  num¬ 
ber  of  wavelengths  in  the  upstream  and  downstream  directions  are  each 
maintained  constant.  The  fj  loop  for  such  a  system  is  shown  in  Figure  2.1. 
The  f2  loop  would  be  similar,  and  the  readout  could  be  the  frequency  differ¬ 
ence  Af  -  £2  -  fp  which  would  be  proportional  to  v.  This  method  may  be 
considered  as  the  cw  analog  of  certain  ultrasonic  pulse  "singaround"  flow¬ 
meters  presently  available  as  commercial  instruments.  ( 15) 

/  2 

Determination  of  v/c 

There  are  a  variety  of  other  methods  of  generating  a  phase  shift  propor¬ 
tional  to  v/  c^.  For  example,  if  one  assumes  that  a  pair  of  opposed  trans¬ 
ducers  could  be  made  identical,  and  that  a  pair  of  identical  circulators 
could  be  perfectly  terminated,  then  uncoded  upstream  and  downstream  cw 
signals  could  be  fed  to  a  differential  amplifier  whose  output  would  be, 
ideally,  proportional  to  the  sine  of  the  phase  angle  between  the  two  waves. 
Offsetting  the  apparent  electronic  simplicity  of  this  method,  however,  is 
the  flow  range  limit  imposed  by  the  lack  of  sensitivity  of  the  sine  function 
for  arguments  near  90°,  as  well  as  the  systematic  errors  and  noise  due  to 
departures  from  the  above  assumptions  on  transducers  and  circulators, 

Still  other  methods  of  separating  upstream  from  downstream  waves  over  a 
common  path  include  rf  bursts  of  fixed  or  variable  frequency,  two  contin¬ 
uous  waves  of  fixed  but  different  frequencies,^  and  cw  alternately  switched 
upstream  and  downstream. 

Some  earlier  investigators  used  symmetrical  dual  path  arrangements  which 
employed  a  single  fixed  frequency  to  drive  the  transducer (s)  radiating  both 
upstream  and  downstream,  and  a  phase  meter  to  sense  the  phase  difference 
between  upstream  and  downstream  receiving  transducers.  ( This  simple 
method  cannot  provide  adequate  absolute  accuracy,  because  of  minor  path 
differences,  but  perhaps  could  be  calibrated  over  a  sufficient  range  of  fuel 
combinations  and  flow  conditions  to  be  useful  in  some  instances. 


Continuous  Wave  Reflection  Densitometer  Method 


Consider  a  main  section  of  buffer  rod  of  impedance  Zg,  with  an  end  plate 
of  impedance  7.j  coupled  to  the  fuel  of  impedance  Z^.  The  specific  acous¬ 
tic  input  impedance  at  the  0  1  interface  is'^) 


Z 


IN 


Z, 


,,  .  ...  .  H.  • 

Z  cos  kd  +  jZ  sm  kd 

C*  _ _ I 

7.  cos  k'cl  !-  jZ,  sin  k"d 

1  u 


(20) 
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Figure  21.  Flow  Velocimeter  Schematic  Showing  One  of 
Two  Loops  for  Maintaining  a  Fixed  Number 
of  Wavelengths  by  Varying  the  Frequency. 


where  the  wave  number  for  the  incident  cw  wave  is  k  =  2tt/X.  and  d  = 
thickness  of  the  end  plate  member.  When  the  plate  is  a  quarter-wave 
thick,  k*d  =  tt/2  and  ZjN  simplifies  to  Z  ^  /  Z^. 


This  permits  one  to  "transform"  the  relatively  low  fluid  impedance  Z2  up 
to  a  value  increased  by  the  factor  Z  ^  /  Z2.  Theory  shows  that  this  can 
substantially  increase  the  sensitivity  of  the  reflection  coefficient  Rq  ^  to 
changes  in  Z^: 


R 


01 


+  Z, 


|Z!2/Z2>  -  Z0 
<2i2/z2)  +  z0 


(21) 


Table  VII  shows  how  Rqj  depends  on  Z2,  the  fluid  impedance,  for  different 
combinations  of  buffer  impedances  Zq  and  Zj.  Theoretically,  this  approach 
could  achieve  a  density  accuracy  of  0.  3%  despite  an  error  in  Rqj  of  2  to  3%. 

For  the  present  application,  however,  we  rejected  this  method  because  of 
complications  to  be  expected  in  a  cw  or  long  rf  burst  approach. 


TABLE  VII.  CW  REFLECTION  COEFFICIENT  RQ1  VERSUS  FUEL 
IMPEDANCE  Z2,  FOR  DIFFERENT  COMBINATIONS 

OF  BUFFER  IMPEDANCES  ZQ  AND  Zl 

Z2 

ZQ=  10 

Z  =  1.  0 

ZQ=  1° 

Z1  =  1.  23 

Zn  =  10 

0 

Z  =  0.  707 

Z  =1.0 

0 

Z  =  0.  316 

Z0  =  6.  67 
Z  =  1.0 

0.  05 

.  333 

.  500 

0 

.  333 

.  500 

.  06 

.  250 

.  430 

091 

.  250 

.  430 

.  07 

.  177 

.  363 

166 

.  177 

.  363 

.  08 

.  Ill 

.  306 

230 

.111 

.  306 

.  09 

.  052 

.  248 

288 

.  052 

.  248 

,  10 

0 

.  200 

333 

0 

.  200 

.  1 1 

047 

,  156 

375 

-.  047 

.  156 

.  12 

091 

.  Ill 

412 

-.  091 

.  Ill 

.  13 

130 

.  074 

444 

-.  130 

.  074 

.  14 

167 

.  034 

-.  473 

167 

.  034 

.  13 

200 

0 

-.500  -.200 

0 
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Indirect  Determination  of  Densit 


Since  p  and  c  are  both  linear  functions  of  fuel  temperature,  it  may  be 
possible  to  determine  p  indirectly  by  measuring  c  and/ or  temperature. 
Figure  22  shows  that  for  fuel  compositions  near  JF-4  this  method  does 
not  appear  sensitive  top.  For  compositions  near  JP-1,  however,  sensi¬ 
tivity  might  be  adequate.  As  a  corollary,  if  the  density  of  virtually  any  un¬ 
known  aviation  fuel  or  fuel  mixture  could  be  independently  measured,  such 
as  with  a  hydrometer,  at  one  or  two  known  temperatures,  the  density1  s 
subsequent  variation  could  be  determined  by  measuring  the  fuel  tempera¬ 
ture.  This  method,  while  undesirable  for  control  purposes,  may  be  ap¬ 
propriate  for  some  diagnostic  studies. 

In  passing,  we  may  note  that  some  present  fuel  control  systems  require 
that  the  fuel  be  known,  such  that  one  of  five  fuel  types  can  be  selected 
manually  on  a  dial  control.  For  such  systems,  where  fuel  type  is  known, 
a  simple  temperature  measurement  would  provide  the  data  for  determining 
density  to  1%  or  better.  For  this  special  case,  an  ultrasonic  flow  veloc- 
imeter  of  the  type  developed  in  this  program  could  provide  an  output  which, 
combined  with  fuel  type  and  temperature  data,  would  measure  M.. 

To  summarize,  and  to  place  the  principal  foregoing  methods  in  perspective 
we  present  Table  VIII.  For  example,  Table  VIII  shows  that  if  one  mea¬ 
sures  flow  velocity  by  a  Doppler  or  phase  method,  then  theoretically  the 
better  choice  for  the  densitometer  is  the  reflection  coefficient  type,  since 
the  outputs  will  be  immune  to,  or  less  sensitive  to,  errors  in  c.  Funda¬ 
mentally,  however,  the  best  of  all  these  combinations  appears  to  be  the 
last.  Here,  both  the  flow  velocity  and  density  mersurcments  are  each 
independent  of  c,  and  naturally,  so  is  their  product,  M. 


TABLE  VIII.  COMPARISON  OF  OUTPUT  PRODUCTS  OF 

SIX  COMBINATIONS  OF  CONTINUOUS  WAVE 
FLOW  VELOCIMETERS  AND  ULTRASONIC 
DENSITOMETERS 

e  n  s  i  to  me  te  r 

Flow  '**> — 

Velocimeter 

Reflection  Coefficient  Resonant  Frequency 
(pc)  (p ) 

Doppler  (v/c) 

Phase  Difference  (v/c^) 
Frequency  Difference  (v) 

M  M/c 

•  /  2 

M/c  M/c 

Me  M 

SOUND  VELOCITY 


TEST  AND  RESULTS 


DOPPLER  FLOW  VELOCIMETER  TESTS 
AT  PANAMETRICS  _ _ 


The  Doppler  effect  was  demonstrated  in  the  test  cell  shown  in  Figures  12 
and  13.  Scattering  of  the  20  MHz  transmitted  beam  was  observed  by 
introducing  one  or  more  small  (diameter  <  1  mm)  gas  bubbles  into  the 
water-filled  test  cell.  As  the  bubbles  rose,  they  scattered  the  incident 
beam  at  90  to  the  receiver  transducer. 

The  Doppler  shifts  were  roughly  estimated  at  2  kHz  from  observations  on 
a  time-domain  oscilloscope.  This  agreed  with  the  bubble  velocity  through 
water.  Since  gas  bubbles  of  ~  5  pm  diameter  have  been  reportedly  de¬ 
tected  by  ultrasonic  waves  even  at  frequencies  below  20  MHz,  it  was  not 
surprising  to  see  large  scattered  returns. (1^) 

To  permit  a  more  quantitative  comparison  between  theory  and  experiment, 
a  spherical  glass  reflector  of  *  0.  1  mm  diameter  was  melted  onto  the  end 
of  a  glass  thread.  When  this  single  reflector  was  moved,  the  correspond¬ 
ing  Doppler  amplitude  was  measured.  Then  the  scattered  power  was  cal¬ 
culated.  As  a  result  of  these  calculations,  it  was  concluded  that  contam¬ 
inant  levels  of  ~1%  of  those  specified  in  MIL-E-5007,  paragraph  3.  4.  1.  3, 
would  be  sufficient  to  yield  detectable  scattered  returns.  It  was  also 
recognized  that  scattering  off  vortices  in  the  fluid  could  turn  out  to  pro¬ 
vide  even  more  scattered  power,  and  might  also  permit  use  of  frequencies 
below  20  MHz.  The  transmission  approach  to  be  described  ultimately 
used  continuous  waves,  whereas  most  previous  transmission  approaches 
used  pulses.  In  our  tests  we  used  both  pulse  and  continuous  wave 
techniques. 

Tests  on  flowing  tap  water,  however,  showed  that  this  fluid  did  not  pro¬ 
vide  detectable  scattering,  using  the  lead  metaniobate  transducers 
available.  Therefore  we  decided  to  build  two  improved  transducers, 
matched,  cut,  ground  and  lapped  to  21  MHz,  1/2  inch  diameter,  of  lith¬ 
ium  iodate.  Even  with  these  improved  transducers,  however,  tests  on 
flowing  water  and  flowing  Avgas -100  did  not  yield  a  sufficiently  high  level 
of  scattered  power  to  assure  meeting  the  accuracy  goal.  At  this  point, 
we  abandoned  the  scattering  approach  in  favor  of  a  transmission  approach. 

TRANSMISSION  FLOW  VELOCIMETER  TESTS 
AT  PANAMETRICS 


Tin:  oscillograms  in  Figures  23  and  24  show  that  in  pulse  transmission 
tests  at  Panametrics  on  water  and  Avgas  - 1  00,  good  signal-to-noise  was 
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obtained,  I'i^uro  2*4  illustrates  an  *»«».•  Ii«»i*  ultrasonic  magnetoatrictive 
wire  experiment  which  simulated  the  separation  of  video  pulses  due  to 
flow.  Subsequently,  the  design  using  rectangular  channels,  and  the  coded 
upstream -downstream  interrogation  method,  was  tested  on  flowing  water. 
The  measuring  cell  was  fabricated  of  2024  aluminum.  Actual  cross- 
sectional  dimensions  of  the  flow  channel  were  0.  502  in.  wide  x  0.  503  in. 
deep.  Cross-sectional  dimensions  of  the  sound  beam  channel  wore  0.251 
in,  wide  x  0,  503  in,  deep,  hi  these  tests  an  Ad- Yu  phase  meter,  model 
422A,  sensed  the  phase  angle,  which  was  proportional  to  v/c>.  The  phase 
meter  was  calibrated  using  a  maximum  water  flow  rate  of  —3.  5  gpm,  or 
—  1750  lb/  hr. 

Figure  2 6  shows  that  the  measured  increase  in  phase  angle  is  approxi¬ 
mately  in  direct  proportion  to  the  volumetric  flow'  rate  of  the  water.  To 
the  extent  that  t ho  temperature  was  constant,  this  graph  shows  that  the 
phase  shift  Is  approximately  proportional  to  M. 

These  tests  also  showed  that  electrical  matching  of  the  transducer  imped¬ 
ances  at  the  circulators  (Figure  15)  was  critical,  and  depended  on  the 
transducer  temperature  and  on  the  acoustic  impedance  of  the  transducer 
hacking  and  that  of  the  liquid.  This  problem  was  partly  overcome  by  ad¬ 
justing  the  matching  network  not  at  zero  flow,  but  at  low  flow. 

These  tests  further  showed  that  new  transducers,  replacing  the  pair  which 
had  boon  originally  electrodod  in  a  rectangular  pattern  to  match  the  inter¬ 
rogation  slot,  could  operate  despite  being  circularly  electrodod.  Thus, 
this  design  simplification  appears  satisfactory. 

When  flow  was  stopped,  the  phase  angle  did  not  return  immediately  to  its 
"zero"  value,  The  delay  appears  tc  be  due  to  curl  or  eddies,  in  the  vi¬ 
cinity  of  the  rectangular  slot,  It  is  expected  that  these  can  be  eliminated 
by  flow  guides  or  acoustic  windows  analogous  to  the  method  of  Figure  6d. 

Following  these  tests,  which  had  been  conducted  at  pressures  below  100 
psig,  the  test  cell,  including  transducers,  was  cyclically  subjected  to 
internal  hydraulic  pressure  up  to  1500  psig.  No  leaks  were  observed.  At 
high  pres  sure ,  minor  fluctuations  in  transmission  were  observed  but  were 
not  considered  to  be  significant. 

DENSITOMETER  TESTS  AT  PANAMETRICS 


Tests  were  made  on  various  r eflectometer -type  probe  materials  (poly¬ 
styrene,  glassy  carbon,  fust'd  silica,  T-40  glass,  AT  quartz  and  titanium), 
on  transducers  (both  hacked  and  unbacked),  and  on  two  electronic  approaches 
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Figure  25.  Signals  Obtained  in  Laboratory-Simulated  Flow, 
Where  Time  Interval  Between  Pulse  Pair  is 
Proportional  to  Distance  That  Coil  is  Away 
From  Center  of  Magnetostrictive  Wire. 
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(logarithmic  comparison  of  rf  bursts  and  linear  comparison  of  the  peaks  of 
spike-excited  pulses).  Densitometer  calibration  tests  used  water,  Avgas- 
100,  acetone,  CS 2  and  CCl^. 

Probe  Materials.  Based  on  the  data  in  Table  VI,  it  appears  that  AT  quartz 
exhibits  the  best  acoustic  properties  among  available  ceramic  materials, 
for  the  probe.  This  material  was  used  in  the  Avco-Lycoming  tests. 

While  the  AT  quartz  densitometer  probe  was  tested  only  near  room  tem¬ 
perature,  we  calculated  the  change  in  probe  impedance  Zj  due  to  changes 
in  density  and  longitudinal  mode  sound  velocity  ever  the  range  -65°  to 
+  160°F  (-54° to+7l°C). 

3 

At  room  temperature,  p  j  -  2.  65  g/  cm  .  The  temperature  coefficient  of 
density  =  -36.  4  x  10“^  per  °C.  Therefore  over  the  full  temperature  range, 
p  ^  changes  by  (71  +  54)  (36.  4  x  10"^)  =  0.  00455  g/ cm^,  or  0.  172%. 


The  temperature  coefficient  of  velocity  Ty  in  a  quartz  crystal  is  given  by 
Mason(^)  in  terms  of  various  elastic  constants  and  the  crystal  orientation 
by  the  following  equation: 

2  2 

T  =  3.  9  +  6.  5  sin  9  cos  </' 
v 

E  ~  E  2  „  2  22 

S11TS11  (cos  ^  cos  ^  +  s*n  ^ 

+  (2s  Ts  +  8  Tsf )  sin^  9  cos^ip(cos^  6  cos^'A  +  sin^tp) 

13  13  44  44 

+  s..Ts..  sin4  9  cos ^\p  -  2s  ^ Ts  sin  9 
3333  1414 


x  sintpcosi^[  3  (cos  9  cos  0  cos  \p  -  sin  0  sin  p)‘ 

2 

-  (sin  <b  cos  9  cos  1 p  4  cos  sirup) _ 2 _ 


£22  22  £22 
s  1  ^  (cos  6  cos  *P  +  sin  1 p)  +  (2a  ^  +  s  )  sin  6  cos  p 

x  (cos^  0  cos^  +  sin^iA)  +  s  ^  ^  sin4  9  cos4  'P 

E  2 

-  2  sin  9  sin^cosip[3(cos0  cos  9  cos  tp  -sin0  sin  \p) 

2 

-  (sin  0  cos  9  cos  i p  +  cos  0  sin  \p )  ] 


22) 
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Retaining  Mason'  s  notation,  we  observe  that  for  an  AT  cut  of  quartz, 

0  =  90°,  &  -  90°,  and  6  =  35°  17'  .  Therefore  most  of  the  terms  in  Eq.  (22) 
drop  out  and  we  find  Tv  =  *2  x  10"^  per  °C,  The  total  change  in  velocity 
is  125  x  2  x  10*6,  or  0.  025%.  As  both  density  and  velocity  coefficients 
are  negative,  the  percentage  impedance  change  is  their  sum:  AZj/Z^  = 


-0.  197%,  or  approximately  +_  0. 
temperature  of  48  F. 


1%  for  -f  112  F  excursions  from  a  mean 


We  measured  the  longitudinal  velocity  in  AT  quartz  as  0.  697  cm/ p  s.  (This 
may  be  compared  with  the  AT  mode  velocity  of  0.  336  cm/ps  corresponding 
to  the  AT  frequency- thickness  constant  of  1.  68  MHz-mm.  )  The  longitudinal 
wave  impedance  Z  j  =  2.  65  x  0.  697  =  1.85  g/  cm**  -p  s  for  AT  qv,  irt  z.  These 
values  appear  in  Table  VI. 


Transducer  Backing.  Backings  were  found  unnecessary,  provided  the 
acoustic  impedances  of  the  probe  and  transducer  were  comparable,  and 
provided  the  bond  between  them  was  sufficiently  thin  and  complete. 


Logarithmic  Comparison  of  Echoes 

Theoretically,  the  apparent  attenuation  (a  L)  increases  as  density  increases, 
where  a  L  =  ln(A/  B),with  A  and  B  defined  in  Figure  17.  Tests  were  run 
with  a  commercial  instrument  which  measures  ln(A/  B),  the  Matec  model 
2470.  These  tests  verified  that  the  model  2470  could  measure  the  ratio  of 
A  to  B  to  +_  0.  0  1  dB,  corresponding  to  B  measurement  accuracy  of  +  0.  1%. 
For  a  buffer  rod  having  Zj<  2  g/cm^-ps,  this  would  give  density  resolu¬ 
tion  of  better  than  +  0.  3%.  That  is,  aluminum,  quartz,  fused  silica  or 
glassy  carbon  exhibit  suitable  impedances. 


Model  2470  stability  tests  showed  drift  and  jitter  *  0.  01  dB  over  110  min¬ 
utes.  But  some  drift,  ~0.  06  dB,  occurred  at  the  beginning  and  end  of  the 
run.  This  drift  probably  could  be  identified  as  to  source,  and  then  reduced 
to  an  acceptable  level.  However,  this  problem,  plus  consideration  of  the 
cost  and  size  of  the  complete  consoles  required  to  operate  the  2470,  includ¬ 
ing  conversion  of  readout  to  density,  led  us  to  concentrate  on  the  alternative 
linear  comparison  of  echoes  approach,  discussed  next. 


Densitometer  Calibration  Tests  in  Various  Liquids 

The  object  of  this  test  was  to  check  our  calculations  of  the  sensitivity  of 
reflection  coefficient  measurements  to  density.  (The  possibility  of  in¬ 
creasing  this  sensitivity  by  measuring  reflection  amplitudes  after  n  re¬ 
verberations  was  also  explored.  )  This  test  was  conducted  using  several 
liquids  of  appropriate  density  and  sound  velocity,  to  simulate  part  of  the 
ranges  experienced  by  Avgas- 100  and  JP-4,  as  listed  in  Table  IX. 
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TABLE  IX.  DENSITY,  SOUND  VELOCITY  AND  ACOUSTIC 
IMPEDANCE  OF  AVGAS- 100  AND  JP-4 


Avgas  - 

100 

JP-4 

Property 

T  - 

-  60°C 

T  = 

+70°C 

T  = 

-60°C 

T  = 

+70°C 

Density  (g/  cm0) 

0.  77 

0.  45 

0.  83 

0.  62 

Sound  Velocity^)  (m/  s) 

1470 

900 

1570 

1070 

Acoustic  Impedance 
(g/cm2-p  s) 

0.  113 

0.  0405 

0.  130 

0. 0664 

m  =  dp  /  dT  (g/  cm^-°C 

0.  00246 

0.  00162 

Liquids  considered  for  simulating  these  ranges  in  room  temperature  tests 
(~25°C)  are  listed  in  Table  X. 


TABLE  X.  PROPERTIES  OF  LIQUIDS  FOR  DENSITOMETER 
CALIBRATION 

Density 

Sound  Velocity 

Acoustic  Impedance 

Liquid 

(g/  cm3) 

(m/  s) 

(g/  cm2-p  s) 

Acetone 

0.  79 

1174 

.  0930 

Carbon  tetrachloride 

1.  595 

926 

.  148 

Kerosene 

.  81 

1324 

.  107 

Water,  distilled 

.  998 

1498 

.  i  50 

Carbon  Disulfide  (25°C) 

1.  26 

1  149 

.  145 

Iso-nctane  (25°C) 

.  687 

- 

- 

Using  the  circuit  of  Figure  20,  and  a  fused  silica  probe,  kA  +  B  was  mea¬ 
sured,  and  density  was  then  calculated  by  the  formula  p  ^  =  (Z  j/ C2>(kA- B)/ 
(kA  +  B)  where  Z  \  =  1 ,  30  g/ cm^ -p,  s .  Results  for  acetone,  water,  CS2  and 
CCI4  are  plotted  in  Figure  27,  covering  densities  from  M).  8  to  1.  6g/ cm’, 
It  is  seen  that  the  data  generally  parallel  the  theoretical  line,  but  errors 
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are  unacceptably  large.  The  errors  appear  largely  due  to  drift  and  noise. 
For  example,  typical  data  for  acetone  were  as  follows:  kA  -  B  -  47  mV, 
and  kA  +  B  =  607  mV.  That  is,  kA  =  327  mV  and  B  =  280  mV.  In  taking 
the  difference  of  these  two  large  numbers,  1  mV  of  noise  produces  a  2% 
error  in  p£.  Since  the  rms  noise  level  for  some  present  broadband 
receivers  is  ~40  p,  V  referred  to  in  the  input,  it  was  expected  that  a  cir¬ 
cuit  such  as  Figure  20  would  not  have  been  subject  to  noise  as  large  as 
1  mV.  Initial  attempts  at  shielding  and  filtering  were  not  sufficient.  More 
work  would  be  required  to  reduce  the  drift  and  noise  to  ~40  p  V,  so  terms 
like  kA  -  B  will  be  repeatable  to  ~0.  1%,  leading  to  errors  in  p  %  of  ~0.  3%. 

PRELIMINARY  TESTS  AT  A VCO- LYCOMING 


Preliminary  transmission  tests  using  video  pulses  were  conducted  at 
Avco-Lycoming  on  February  22,  1972.  The  fluid  was  Stoddard  solvent, 
MIL-F7024A,  Type  2.  Mass  flow  rate  was  varied  from  zero  up  to  2800 
lb/  hr  at  35  psig.  The  ultrasonic  test  cell  is  shown  in  Figure  28  (except 
for  Aeroquip  pipe/ AN  adapters,  2021-8-18).  Because  of  the  long  ultra¬ 
sonic  transmission  path,  12  in.  ,  signal  identification  was  poor  at  low  flow. 
But  at  turbulent  high  flow,  pulse  separation  was  clearer.  Pulse  separa¬ 
tions  were  similar  to  those  of  Figure  23.  Results  appeared  to  be  similar 
to  those  obtained  earlier  at  Panametrics  using  Avgas- 100  (Figure  29), 
circulated  in  a  closed  loop  (Figure  30). 

To  determine  the  effect  of  engine  noise  on  the  ultrasonic  transmission 
measurement,  the  ultrasonic  unit,  still  filled  with  the  calibrating  fluid, 
was  strapped  with  wire  to  a  T43  LI  3  full  engine.  The  unit  was  located  next 
to  an  accelerometer,  directly  over  the  fuel  control  part  of  the  engine,  to 
subject  it  to  realistic  vibration  levels.  The  engine  was  then  run  for  about 
5  minutes  at  speeds  up  to  97%,  Rig  vibration  was  measured  as  0.  001  in. 
(double  amplitude,  filtered  below  600  Hz).  Using  a  6  MHz  portable  oscil¬ 
loscope,  the  noise  was  observed  to  be  <  50  mV  peak-peak,  at  ~100  kHz, 

The  no-flow  10  MHz  signal  was  200  mV  peak- peak,  for  a  pulse-echo  total 
path  of  24  in. 

TESTS  ON  JP-4,  JP-5  AND  A  CONTAMINATED  FLUID 

The  main  purpose  of  these  tests,  conducted  at  Avco-Lycoming  between 
June  and  August  1972,  was  to  determine  the  accuracy,  reproducibility  and 
linearity  of  the  ultrasonic  mass  flowmeter  system,  and  to  demonstrate 
operation  on  uncontaminated  fuel  as  well  as  on  fuel  contaminated  in  accord¬ 
ance  with  MIL-E-5007C,  paragraph  3,  4.  1.  3,  The  response  time  was  not 
determined  in  these  tests,  but  was  calculated  from  the  system  bandwidth 
to  be  20  ms. 
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RANSIT  TIME  DIFFERENCE ,  ns 


6  6  gpm  at  ZOOOrpm 
12  Ogpm  at  3400rpm 


Figure  30.  Flow  Loop  for  Testing  Scattering  and 
Transmission  Methods  of  Measuring 
Flow  Velocity, 


Description  of  Test  Stand.  The  flowmeter  calibrating  stand  was  manufac¬ 
tured  by  Cox'  s  Commercial  Research  Division,  the  George  L.  Nankervis 
Company,  Detroit.  It  is  denoted  Part  No.  T-303T,  Serial  No.  6602-A, 
and  includes  readouts  such  as  a  timer,  weigh  switch,  hydrometer,  vis¬ 
cometer,  fuel  temperature,  fuel  pressure,  and  was  used  with  two  turbine- 
type  flow  sensors  to  cover  the  flow  ranges  of  interest.  Figures  31  and  32 
illustrate  the  calibration  test  arrangement. 

The  following  results  were  obtained  in  JP-4and  JP-5  prior  to  contaminant¬ 
testing  of  the  ultrasonic  measuring  cell. 

Calibration  Results  on  JP-4.  Calibration  data  for  JP-4  is  given  in  Table 
XI  for  the  following  conditions:  Date:  June  12,  1972,  Cox  turbine  flow¬ 
meter  Serial  No.  8463,  size  AN8-6,  range  188.  4  to  2062.  8  lb/hr.  Tem¬ 
perature:  76°F.  Fuel  specific  gravity:  0.  760.  Kinematic  viscosity:  1.  04 
centistokes.  The  ultrasonic  transit  time  across  a  1-in.  path  perpendicular 
to  the  flow  was  20.  85  (is,  measured  with  the  5220  gage.  Thus  the  sound 
speed  was  0.  480  in.  /  p  s  or  1220  in/  s,  and  the  fuel1  s  characteristic  acous¬ 
tic  impedance  was  =  0.  0926  g/  cm“-p  s,  or  about  62%  of  that  of  water 
at  the  same  temperature. 


TABLE  XI.  FLOWMETER  CALIBRATION 
FOR  JP-4 

TEST 

DATA 

Cox 

Cal.  Pts. 

(Hz) 

Avg.  Flow 

Cal.  Wt.  Time 

(lb)  (s) 

Computed  Ultrasonic  Readings  of 

Mass  Flow  Rate  Phase  Angle  Density  M 
(lb/hr)  (Deg)  (Relative  Units  on  DVM] 

1320 

1 1 

19.  197 

2062.  8 

25 

862 

2409 

3000 

1200 

10 

19.  183 

1876.  7 

24 

817 

2411 

2800 

1080 

9 

19.  2  i  3 

1686.  4 

22 

750 

2415 

2550 

900 

7.  5 

19.  197 

1406.  5 

15.  5 

510 

2400 

1670 

720 

6 

19.  147 

1128.  1 

13 

4  30 

2400 

1200 

600 

5 

19.  120 

941.  4 

7 

250 

2407 

499 

480 

4 

- 

- 

- 

- 

- 

- 

360 

3 

19.  027 

567.  6 

2.  5 

100 

2409 

181 

300 

2.  5 

- 

- 

- 

- 

• 

- 

180 

1.  5 

19.  133 

282.  2 

- 

- 

- 

- 

120 

1 

19.  1  10 

188,  4 

- 

- 

- 

w 

Figure  1|.  Block  l.)iaj>  ram  of  Ca  1  il>  rat  ion  list  A  nMmji'Hn-nli 


Figure  32.  View  of  Ultrasonic  Equipment 
and  Calibration  Test  Stand. 
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The  ultrasonically- determined  M  is  plotted  against  the  standard  computed 
mass  flow  rate  in  Figure  33.  It  is  seen  that  the  results  appear  to  be 
nearly  linear,  with  deviations  from  a  straight  line  fit  of  about  100  lb/hr, 
or  +_  5%  of  range.  More  data  are  needed,  however,  for  better  curve 
definition. 

Besides  the  JP-4  data  tabulated  above,  additional  measurements  were  ob¬ 
tained  of  phase  angle  versus  Cox  M  calibration  points.  These  data  are 
plotted  in  Figure  34. 

Results  are  similar,  but  not  identical.  Drift  or  scatter  may  be  mainly 
accounted  for  by  the  phase  meter.  In  future  tests,  phase  accuracy  of  0.  3 
would  eliminate  this  possible  source  of  error. 

Calibration  Results  on  JP-5.  Calibration  data  for  JP-5  is  given  in  Table 
XII  for  the  following  conditions:  Date:  June  13,  1972.  Cox  turbine  flow¬ 
meters,  Serial  Nos.  8463  (as  before)  and  5903,  size  AN10,  range  817.  1  to 
4957.  9  lb/hr.  Temperature:  74°F.  Pressure:  22  psig.  Fuel  specific 
gravity:  0.  820.  Kinematic  viscosity:  2  centistokes.  The  ultrasonic 
transit  time  across  a  1-in.  path  perpendicular  to  the  flow  was  19-24  p  s, 
measured  with  the  5220  gage.  Thus  the  sound  speed  was  0.  052  in.  /ps  or 
1320  m/s,  and  the  fuel's  acoustic  impedance  was  2.^  -  0.  108  g/cm^-ps, 
or  about  72%  of  that  of  water  at  the  same  temperature. 

The  ultrasonically-determined  M  is  plotted  against  the  standard  com¬ 
puted  mass  flow  rate  in  Figure  35.  Significantly,  for  3  of  4  runs  in 
Figure  35,  the  data  are  within  +.  4%  of  one  straight  line.  Disregarding 
the  anomalous  fourth  run,  these  data  are  quite  encouraging,  in  that  they 
demonstrate  the  basic  linearity  and  reproducibility  of  the  M  system.  In 
Figure  36,  where  Ivl  ranged  to  5000  lb/hr,  beyond  the  range  of  our  multi¬ 
plier,  we  show  that  the  measured  phase  angle  may  be  approximated  as  a 
piecewise  linear  function.  Departures  from  linearity  are  apparently  due 
to  eddies.  . 

Contamination  Tests.  The  stand  used  in  these  tests  was  manufactured  by 
Engineered  Devices  of  Agawam,  Massachusetts.  Tests  were  run  at  ~1900 
lb/hr,  as  measured  by  the  Cox  Series  12  rig  rotameter.  The  fluid,  con¬ 
taminated  essentially  in  conformance  with  MIL.-E- 5007C,  paragraph 
3.  4.  1.  3,  was  introduced  without  filtering  at  ~35  psig  and  ~75°F,  with  at 
least  one  brief  temperature  rise  of  ~105  F.  Tests  began  on  June  14  and 
ended  on  August  24,  1972. 

The  MIL-E-5007C  types  and  distributions  of  contaminants  are  given  in 
Table  XIII.  (Instrumentation  specification  MIL-F- 86l  5B  is  the  same  except 
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Figure  Ullrusonicnlly  DHcrmincri  fvl  Cor  JP-4. 

(>') 


Lr 


6  8  iOO  0  iZ  J4 

JP-4  CALIBRATED  MASS  FLOW  RATE.  LB  /HR 


ULTRASONIC  PHASE  SHIFT.  RELATIVE  UHITS 


ULTRASONIC  ALLY  DETERMINED  MASS  FLOW  RATE. RELATIVE  OMITS 


ULTRASONIC  PHASE  ANGLE, DES 


Figure  36,  Ultra^onically  Determined  Phase  Angle 
Vers  js  Xvl  for  JP-5. 
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TABLE  XII.  FLOWMETER  CALIBRATION  TEST  DATA 
FOR  JP-5 


Cox 

Cal.  Pts. 

(Hz) 

Cal.  Wt. 

(lb) 

Avg.  Flow 

T  ime 

(s) 

Computed 

Mass  Flow  Rate 

(lb/ hr) 

Ultrasonic  Readings  of 

Phase  Angle  Density  M 

(Deg)  (Relative Units  on  DVM) 

1320 

1 1 

17.  710 

2236 

84 

2700 

2175 

12000 

1200 

10 

17.  690 

2035 

66.  5 

2152 

2134 

7650 

1080 

9 

17.  680 

1833 

56 

1830 

2131 

8350 

900 

7.  5 

17.  680 

1527 

41 

1233 

2133 

5400 

720 

6 

17. 687 

1221 

28 

925 

2139 

4150 

600 

5 

17.  700 

1017 

19.  5 

650 

2127 

2850 

480 

4 

17.  693 

814 

1 1 

388 

2141 

1630 

360 

3 

17.  740 

609 

4 

180 

2150 

750 

300 

2.  5 

17.  713 

508 

3 

214 

2148 

548 

180 

1.  5 

17.  800 

303 

- 

- 

- 

- 

120 

1 

17.  710 

203 

- 

- 

- 

- 

(A  bove : 

Cox  S/N  8463) 

(Below: 

Cox  S/N  5903) 

1220 

25 

17.  853 

5041 

145 

4580 

2172 

13061 

1200 

25 

18.  153 

4958 

142 

4510 

2174 

13061 

1 100 

25 

19.  827 

4539 

124 

3961 

2179 

13062 

1000 

25 

21.  817 

4125 

103 

3251 

2180 

13061 

900 

25 

24.  270 

3708 

84 

2660 

2177 

1 1400 

800 

25 

27.  307 

3296 

66 

2127 

2198 

8893 

700 

25 

31.  250 

2880 

49 

1566 

2198 

6200 

600 

25 

36.  490 

2466 

38 

1198 

2183 

4387 

500 

10 

17.  550 

2051 

27 

851 

2190 

2882 

400 

10 

21.  98  3 

1638 

16 

499 

2181 

1109 

300 

10 

29.  330 

1227 

2.  5 

85 

2189 

395 

200 

10 

44.  060 

817 

- 

- 

- 
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the  concentrations  of  the  first  four  types  are  increased  by  2.  5  by  requir¬ 
ing  the  same  quantities  per  400  gallons  instead  of  per  1000  gallons.  ) 


TABLE  XIII.  FUEL  ENDURANCE  TEST 

CONTAMINANT 

Contaminant 

Particle  Size 

Quantity 

Iron  oxide 

0-5  microns 

5-10  microns 

28.  5  gr>  '1000  gallons 
1.  5  gm/  1000  gallons 

Sharp  silica  sand 

150-300  microns 

300-420  microns 

1.  0  gm/  1000  gallons 
1.  0  gm/  1000  gallons 

Prepared  dirt  con¬ 
forming  to  AC  Spark 
Plug  Co.  Part  No. 
1543637  (coarse 
Arizona  road  dust) 

Mixture  as  follows: 

0-5  microns  (12%) 

5-i0  microns  (12%) 
10-20  microns  (14%) 
20-40  microns  (23%) 
40-80  microns  (30%) 
80-200  microns  (9%) 

8.  0  gm/  1000  gallons 

Cotton  linters 

Grade  6,  staple  below  7, 
second  cut  linters  (U.  S. 
Department  of  Agri¬ 
culture  Grading 
Standards). 

0.  1  gm/1000  gallons 

Crude  naphthenic  acid 

0.  03%  by  volume 

Salt  water  solution  shall 
contain  4  parts  NaCl 
to  96  parts  H^O  by  wt. 

0.  01%  entrained 

The  tests  were  originally  scheduled  to  continue  on  a  one-shift  basis  for 
100  hours.  However,  loss  of  signal  at  25  hours  caused  an  interruption. 
Examination  of  the  test  cell  and  transducers  at  Avco- Lycoming  (Figures 
37,  38)  showed  some  discoloration  and  corrosion,  and  some  minor  slurry 
accumulation,  but  not  enough  to  account  for  loss  of  signal.  The  complete 
system  was  returned  to  Panametrics,  where  it  was  found  that  the  5  MHz 
frequency  had  drifted  outside  the  narrow  range  of  the  receiver  circuitry. 
Retuning,  and  readjustment  of  the  code  delay,  restored  the  signals. 


Prior  to  returning  the  system  to  A  v'co-  Lycoming ,  it  was  endurance-tested 
for  100  hours  on  unfiltered  tap  water  at  ~800  1b/hr,  plus  overnight  and 
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weekend  soaks  at  no-flcw,  with  periodic  calibration  checks  against  an 
accumulated  volume  of  water.  The  flow  velocimeter  appeared  quite  re¬ 
producible  . 

The  system  was  returned  to  Avco-Lycoming  and  contamination  testing  was 
resumed  on  August  10,  continuing  until  August  24,  1972.  At  this  termina¬ 
tion,  exposure  totaled  104  hours  at  1840  to  1940  lb/hr,  plus  16  hours  of 
filtered  contaminant  at  ~1900  lb/  hr,  plus  second-  and  third- shift  and  week¬ 
end  soaks  at  no-flow  totaling  428  hours.  Signal  levels  decreased  in  these 
tests  from  a  maximum  of  2V  to  a  minimum  of  0.  06V.  This  decrease  how¬ 
ever  appears  to  be  mainly  a  result  of  amplifier  malfunction,  not  transducer 
degrading.  Since  phase  meters  are  now  available  which  maintain  their  ac¬ 
curacy  for  signals  above  0.  01V,  this  decrease  of  signal  amplitude  is  not 
considered  significant. 

The  ultrasonic  measuring  cell  was  opened  for  inspection  prior  to  recalibra¬ 
tion.  It  was  observed  that:  the  transducer  vertical  faces  were  partly 
covered  with  a  dark  red  residue;  the  densitometer  probe'  s  horizontal  face 
showed  one  red  speck;  some  dark  red  slurry  had  accumulated  at  the  bot¬ 
tom  of  the  transducer  cavities;  and  some  corrosion  and  deposits  were 
evident  at  mating  surfaces.  The  mating  surfaces  were  wiped  clean  to  ease 
reassembly,  and  then  the  cell  was  tapped  together  again.  No  residue  was 
removed  from  the  flow  channel  nor  from  the  ultrasonic  measuring  channel. 

Next,  the  system  was  transferred  back  to  the  Cox  calibration  test  stand 
(Figures  31,  32)  and  then  recalibrated  on  JP-4  at  77.  5°F  and  at  92°F. 
Figure  39  shows  the  measured  phase  difference  versus  M,  obtained  using 
the  flow  velocimeter  portion  of  Figure  15.  At  77.  5  F  we  also  measured 
the  phase  difference  versus  M  by  manually  switching  the  uncoded  5  MHz 
wave  alternately  upstream  and  downstream.  This  method  is  valid  to  the  ex¬ 
tent  that  the  average  c  remains  constant  for  both  directions.  Note  that  the 
use  of  coded  waves  continuously  transmitted  both  ways  avoids  this  limita¬ 
tion,  and  is  obviously  to  be  preferred. 

Recalibration  Results  on  JP-4,  Recalibration  data  for  JP-4  is  given  in 
Table  XIV  for  the  following  conditions:  Test  Date:  August  24,  1972.  Cox 
turbine  flowmeter  Serial  No.  8932,  size  AN8,  range  214  to  3561  lb/hr, 
calibrated  March  20,  1972.  Test  temperature:  77.  5  F,  Fuel  specific 
gravity:  0.  751. 

The  essential  corresponding  data  for  JP-4  at  92°F  is  given  in  Table  XV. 

Fuel  specific  gravity  was  0,  746.  Coded  5  MHz  waves  were  continuously 
transmitted  both  ways. 
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PHASE  ANGLE,  DEGREES 


Figure  39.  Recalibration  Curves  for  JP-4  at  Different 
T  emperatures. 
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TABLE  XIV.  RECALIBRATION  DATA  FOR  JP-4 
AT  77.  5°F 


Cox 

Cal,  Pts. 

(Hz) 

Cal.  Wt. 

(lb) 

Avg,  Flow 

T  ime 

(s) 

Computed 

Mas  s  Flow  Rate 

(lb/hr) 

Phase  Angle 

5  MHz  Coded  CW  Alter 

Run  l  Run  2  Run  3  Run  4 

(degrees) 

Dir. 

Run  5 

90 

2.  25 

37.  857 

214.  0 

_ 

_ 

120 

3.  0 

37.  887 

285.  1 

- 

- 

- 

- 

180 

4.  5 

37.  847 

428.  0 

- 

- 

- 

- 

- 

240 

6.  0 

37.  897 

570.  0 

- 

- 

- 

- 

300 

7.  5 

37.  957 

711.  3 

- 

- 

- 

10 

- 

360 

9.  C 

37.  937 

854.  0 

- 

- 

- 

13 

14 

480 

12.  0 

37.  860 

1  141.  0 

4 

5 

5 

25 

- 

600 

15.  0 

37.  947 

1423.  0 

8 

6 

6 

38 

40 

720 

18.  0 

37.  940 

1708.  0 

17 

13 

12 

45 

- 

900 

22.  5 

37.  923 

2135.  9 

30 

26 

25 

64 

59 

1080 

27.  0 

37.  923 

2563.  1 

42 

42 

41 

73 

78 

1200 

30.  0 

37.  003 

2849.  4 

53 

53 

50 

- 

80 

1500 

37.  5 

37.  910 

3561.  1 

75 

75 

73 

- 
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TABLE  XV.  RECALIBRATION  DATA 
AT  92°F 

FOR  JP-4 

Computed 

Phase  Angle 

Mass  Flow  Rate 

Run  6 

Run  7 

Run  8 

(lb/  hr) 

( degrees) 

210 

- 

280 

- 

420 

- 

560 

- 

698 

- 

838 

- 

. 

1120 

1  1 

1  1 

1  1 

1397 

13 

16 

17 

1677 

23 

25 

30 

2096 

40 

41 

52 

25  16 

56 

58 

72 

2797 

68 

68 

84 

3504 

14 

04 

1  10 
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In  Figure  39,  the  three  runs  of  v/c  or  phase  data  for  77.  5°F  show  minor 
hysteresis  or  drift,  which  amounts  to  +  3%  of  the  range.  For  92°F,  runs 
6  and  7  show  even  less  spread,  generally  less  than  +  2%.  Run  8  appears 
anomalous,  and  is  probably  due  to  a  fuel  temperature  increase.  Diffi¬ 
culty  was  encountered  in  controlling  fuel  temperature  during  these  tests. 
In  summary,  these  recalibration  data  points  show  that  the  cell  and  trans¬ 
ducers  were  able  to  function  after  prolonged  exposure  to  a  contaminated 
fluid. 

Figure  40  demonstrates  that  if  one  does  not  transmit  simultaneously  in 
both  directions,  scatter  of  the  data  increases,  since  small  variations  in 
c  influence  the  v/c^  readout.  These  results  involved  the  manual  switch¬ 
ing  or  alternating  of  the  propagation  direction,  with  phase  being  averaged 
over  several  seconds  for  each  direction. 


Figure  40.  Phase  Angle  Difference  for  Alternate 
Propagation  Directions  for  JP-4. 
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SUMMARY  OF  RESULTS 


The  principal  results  of  the  theoretical,  experimental,  and  test  portions  of 
this  program  may  be  summarized  as  follows: 

•  The  ultrasonic  mass  flowmeter  measuring  cell  has  been  demonstrated 
to  operate  without  clogging,  both  during  and  after  exposure  to  a  fluid 
contaminated  essentially  in  conformance  to  MIL-E-5007C.  Exposure 
included  shifts  of  approximately  8  hours  each,  totaling  104  hours  at 
approximately  1900  lb/hr,  plus  overnight  and  weekend  soaking  of  428 
hours  at  no-flow.  Endurance  appears  to  be  far  in  excess  of  the  two 

2.  5-hour  tests  required  for  instruments  in  MIL-F-3615. 

•  A  pseudo- random-noise  phase  coding  technique  has  been  demonstrated 
to  provide,  for  the  first  time,  separation  of  two  continuous  waves  of 
the  same  frequency  transmitted  upstream  and  downstream  over  a 
common  path.  This  technique  provided  a  flow  velocimeter  response 
time  of  20  ms,  a  factor  of  25  faster  than  the  contractual  objective  of 

0.  5  s,  and  a  range  up  to  5000  lb/hr.  It  was  found  that  undesirable 
circuit  unbalances  were  produced  by  fluctuations  in  fuel  'emperaturo 
of  a  few  degrees,  or  by  changes  in  the  acoustic  radiation  impedance 
due  to  changes  in  impedance  of  the  fuel  or  transducer  backing.  To 
overcome  this  problem,  it  will  be  desirable  to  choose  a  transducer 
with  minimal  variation  in  properties  due  to  temperature,  and  to 
isolate  the  transducer  from  the  fuel,  using,  for  example,  an  attenu¬ 
ating  acoustic  window  as  a  pad. 

•  Immunity  to  uncertainties  in  Reynolds  number  or  uni  e rtainties  in 
unidirectional  flow  profiles  appears  achievable  using  the  square  (low 
channel  wherein  1  00"/>,  of  the  fuel  cross  section  is  interrogated  by  a 
rectangular  plane  wave  beam  of  ultrasound.  This  is  the  basis  lor 
accurate,  linear  response  to  flow  velocity,  independent  of  luel  com¬ 
position,  temperature,  density  or  viscosity.  In  the  present  measur¬ 
ing  cell,  the  observed  departures  from  linearity  appear  to  be  due  to 
nonideal  inlet  conditions,  probably  correctable  with  a  flow  straightener 
and  to  fuel  rotation  at  the  measuring  region,  probably  correctable 
with  flush-mounted  acoustic  windows,  wedges,  screens  or  other  means 

•  Based  on  the  encouraging  first  results  obtained  with  coded  single  fixed 
frequency  waves,  t  appears  that  the  use  of  two  variable  1  requeue ies 
should  work  too,  to  measure  v  independent  of  c.  This  is  the  basis  tor 
in  improved  ultrasonic  mass  flowmeter,  where  both  v  and  are  each 
measured  independent  of  c.  The  product  of  their  outputs  is  thus  pro¬ 
portional  to  M, 


•  Since  signal  levels  in  some  cases  were  as  small  an  20  mV  rms,  Anti 
whereas  the  accuracy  of  the  Ad-Yu  phase  meter  available  fur  the  tests 
In  spec  Hied  fur  signals  above  300  mV  nun,  il  In  understood  that 
much  ui  the  xyitcm'it  present  inability  it*  measure  Kf  with  sufficient 
accuracy,  especially  below  ^  800  lb/hr,  corresponding  to  relatively 
small  phase  angles,  In  attributable  it*  operating  beyond  thin  meter*  s 
limitations,  Kurt  bur,  (In*  (act  (hat  tin*  extrapolated  results  do  nut 
pANN  through  tin*  origin,  in  Figure*  26,  33,  34,  38,  30  untl  40,  in 
utt ributed  it*  tin*  manual  "aero"  adjustment  of  the  phase  meter  at  a 
non- aero  flow,  This  adjustment  wan  required  in  order  to  match  the 
present  circuits  to  tin*  t  ransduce i* ,  whose  electrical  impedance  de¬ 
pend*  on  the  acouNtic  impedance  of  tin*  adjacent  media,  and  on  it* 
own  temperature.  It  In  recognised,  however,  that  phase  muter  re¬ 
quirements  such  as  automatic  0  adjust,  0, 3U  accuracy  on  10  mV  sin- 
unis,  faster  warmup  aim  10  ms  response  can  now  be  met  by  a  new 
instrument  available  commercially,  Also,  uni*  of  a  different  tem¬ 
perature  -  immune  pieaoeluct ric  crystal,  instead  of  a  ceramic  trans¬ 
ducer,  and  acoustical  isolation  or  padding  of  this  crystal  from  (he 
variable  load  presented  by  different  fuels,  is  expected  to  avoid  the 
other  observed  deficiencies,  Accordingly,  range,  accuracy,  response 
time  and  other  goals  ut'e  considered  achievable. 


*  For  the  Ad-Yu  422A,  the  relative  accuracy  la  1/4°  plus  a  high  fre¬ 
quency  error  of  0.  lu  per  MHss,  The  absolute  error  is  +_  1°  n  ±  2%, 
whichever  is  greater.  In  addition,  there  is  a  high  frequency  error 
of  0,  1°  per  Mils  or  0.  1%  of  the  reading  per  MHa,  whichever  is 
greater.  For  tests  at  5  MHn,  this  amounts  to  _+  0,  5%,  for  signals 
between  300  mV  and  20  V  rma.  A  further  error  source  especially 
at  low  flow  rates  (small  phase  angles)  is  drift.  For  example,  after 
one  hour  warm  up,  the  meter  reading  will  drift  less  than  lu  in  ten 
hours.  This  drift  must  be  reset  manually  with  the  panel  zero  set 
control. 
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CONCLUSIONS 


All.' r  examining  the  program1  h  results,  live  major  conclusions  may  lu> 
drawn! 

•  A  now  ultrasonic  transmission  flow  velocimeter  technique  for  con- 
ttnuously  measuring  v/c*-  has  boon  demonstrated.  Thin  technique 
»s  believed  lo  inherently  offor  the  best  accuracy  consistent  wltl\  fast 
respunse,  20  mti,  of  any  ultrasonic  approach  demonstrated  so  far. 
Operation  up  lo  *>000  lb/ hr  was  demonstrated, 

•  A  now  ultrasonic  reflection  densitometer  technique  for  measuring  pc 
has  boon  demonstrated  on  liquids  ranging  in  density  from  about  0.  7 
to  1 .  6  g/  cm 

•  A  now  method  of  interrogating  100%  of  the  flowing  liquid1  s  square 
cross  section  using  a  rectangularly  collimated  plane  wave,  which 
potentially  offers  linearity  and  accuracy  essentially  independent  of 
Reynolds  number  or  unidirectional  flow  profile,  lias  been  demon* 
strated  on  water,  Avgas- 100,  JP-4,  JP-ii  and  a  contaminated  solvent. 

•  Using  a  cell  design  containing  nonintrusive  trans»...eer s  and  no  rotat¬ 
ing  parts,  endurance  of  over  100  hours  has  been  demonstrated  using 

a  test  fluid  flowing  at  ~1‘>00  lb/  hr  and  contaminated  essentially  in  con¬ 
formance  with  military  specifications  for  flow  measuring  instruments. 

•  A  new  ultrasonic  muss  flowmeter  has  been  built,  consisting  of  a  flow 
velocimeter  (v/c^),  a  densitometer  (pc)  and  a  time  intervalomoter 
(T  ■-  l/c),  The  combination  of  these  outputs  has  been  shown  to  yield 
an  output  voltage  approximately  proportional  to  the  mass  flow  rate  Ivl 
over  specified  ranges,  Departures  from  linearity  are  presently 
mainly  attributed  to  stationary  eddies  which  are  partly  in  the  measur¬ 
ing  path.  It  is  expected  that  they  can  be  virtually  eliminated  using 
acoustic  windows,  wodges,  screens  or  other  means. 

Based  on  the  theoretical,  experimental  and  test  results,  it  is  concluded 
that  most  of  the  observed  inaccuracy  in  v/  c^  (about  _f  3%  at  high  flow)  and 
in  p  c  (about  +  5%)  and  therefore  nearly  _+  10%  in  Ivi  at  2000  lb/ hr  can 
be  attributed  to  factors  which  can  be  eliminated.  For  example  the  in¬ 
accuracy  of  the  phase  meter  available  in  this  program  to  measure  v/ 
was  +  2°,  Phase  meters  now  available  can  reduce  this  source  of  inaccu¬ 
racy  to +.0,1°.  Further,  in  the  pc  cast',  the  observed  drift  and/or  noise 
of  several  mV  produces  an  error  in  p  of  2%  per  mV,  yet  drift  and  noise 
levels  of  0,  1  mV  or  less  are  considered  achievable.  Thus,  p  c  impr ecision 
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Of  i' f f of 8  no  larger  thun  0,  2%  appear  achievable.  ll  is  therefore  esti¬ 
mated  that  a  (v/e^)  (pc)  (l/T)  system  could  achieve  an  Ml  accuracy  of 
belief  than  1%  at  2000  lb/hr,  and  as  a  corollary,  a  pv  system  indepen¬ 
dent  of  c  should  provide  even  better  accuracy,  since  c  errors  are  elim¬ 
inated,  and  other  electronic  computation  error  sources  would  be  reduced 
in  number, 

Several  other  conclusions  l'rom  this  work  include: 

1.  Any  type  of  ultrasonic  measurement  in  the  t'nr  field  is  subject  to 
serious  error  due  to  wavefront  curvature  under  laminar  flow, 

2>  In  the  near  field,  time  of  flight  measurements  using  pulses  woulc, 
require  subnanosecond  accuracy,  and  fast  response,  beyond  the 
limits  of  present  instruments, 

3,  Doppler  measurements  could  properly  Weight  the  profile  if  inson- 
ification  were  uniform  and  if  scatter ers  were  uniformly  distributed. 
Alternatively,  range  gating  of  waves  obliquely  incident  upon  axi- 
symmetric  flow  could  weight  the  profile,  Uniformity  of  scatterers 
could  be  enhanced  using  a  Kenics  static  mixer  containing  at  least 

7  sections,  which  is  probably  an  acceptable  number  based  on  length 
and  pressure  drop  for  a  1/ 2-inch-diameter  pipe.  The  dual-frequency 
Doppler  system  would  sense  v/  c  in  the  scattering  -/.one,  and  would  be 
relatively  immune  to  v/  c  outside  that  zone.  If  a  static  mixer  is  used, 
sensing  should  be  accomplished  within  ~1  to  3  diameters  of  the  mixer 
exit.  High  Doppler  S/N  obtained  on  a  standard  cross  section  (glass 
sphere)  and  on  bubbles  implied  that  useful  Doppler  measurements  on 
fuels  contaminated  above  ~1%  of  MIL-E-5007C  may  be  achievable, 
but  accuracy  and  response  time  may  be  sacrificed,  especially  on 
uncontaminated  fuel. 

4,  If  a  clamp-on  flowmeter  is  to  be  installed  on  a  standard  all-metal 
straight  tube  or  pipe,  conical  waves  can  provide  a  measure  l'  vj 
independent  of  Reynolds  number  or  flow  profile.  Difficulties  anti¬ 
cipated  in  measuring  c  and  p  for  this  geometry  suggest  that  a  two- 
frequency  technique  be  employed,  to  yield  v^  (not  v^/  c  or  v(j/  c^).  p 
data  may  be  obtained  nonacoustically  (e.  g,  ,  from  fuel  type  and  tem¬ 
perature  inputs)  if  Ml  is  the  desired  final  output. 

Analysis  of  v  versus  a  dons ity- temperature  function,  taken  with  the 
profile  equations,  shows  that  if  p  and  T  are  known  to  1%,  v  and 
hence  Re  can  be  determined  to  20%  or  belter,  for  which  the  uncer¬ 
tainty  in  K  =  va/vc|  be  <  0.  1  % ,  if  the  profile  were  tully 
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developed.  On  the  other  hand,  determination  of  K  from  measure- 
ments  of  v(j  alonr  is  generally  not  possible  to  better  than  1%,  since 
in  the  present  problem  the  viscosity  .and  Re  ranges  span  over  a 
factor  of  10  for  different  (unknown)  fuels  even  <il  known  temperatures. 

h.  The  Doppler  method,  while  not  used  in  the  final  flowmeter,  might 
nevertheless  be  useful  to  detect  contaminants  well  below  1%  of  M1L- 
E-8007C,  or  to  measure  v/c  (average  and/  or  profile)  in  fluids  in 
which  seutterers  such  as  bubbles  or  microballoons  could  be 
introduced. 

The  above  conclusions  may  be  summarized  by  observing  that  a  refinement 
of  the  new  technology  demonstrated  in  this  program  should  result  in  an 
ultrasonic  mass  flowmeter  capable  oi  meeting  current  USAAMRDL.  re¬ 
quirements  for  diagnostics  and  control. 
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RECOMMENDATIONS 


Baaed  on  the  present  work,  it  appears  that  the  present  ultrasonic  mass 
flowmeter  system  can  be  refined.  This  new  system,  without  deleting  any 
of  the  essential  requirements  of  the  present  program,  could  additionally 
provide  a  mass  flow  reading  which  is  substantially  a  linear  function  of  M 
with  a  flow  velocity  response  time  of  20  ms,  and  in  one  arrangement  may 
provide  v  and  p  measurements  each  independent  of  c.  To  achieve  this 
improved  performance,  we  recommend  the  following  approach: 

•  Optimize  the  present  ultrasonic  mass  flowmeter  to  combine  v/c~, 
pc  and  T  with  cell  improvements  such  that  the  output  is  essen¬ 
tially  in  direct  proportion  to  M. 

•  Determine  the  limits  of  a  commercially  available  resonant  vane 
densitometer,  particularly  with  respect  to  contaminants  and 
vibration. 

•  Assuming  the  resonant  densitometer  (p  )  device  operates  satis¬ 
factorily,  add  to  the  v/c^  equipment  a  two-frequency  capability 
such  that  v  can  be  determined  independent  of  c. 

•  Optimize  the  combination  of  p  and  v  to  provide  M,  where  p  and 
v  are  each  independent  of  c. 

•  Conduct  preliminary  calibration  and  endurance  tests  using  water 
and  Avgas  flowing  at  rates  of  up  to  ~2000  lb/hr,  and  demonstrate 
densitometer  and  transducer  operation  from  -65°  to  +  160°F,  at 
least  under  no-flow  conditions. 

•  Conduct  calibration  tests  and  contaminant  tests  as  per  the  latest 
applicable  military  specification  on  instrument  endurance.  The 
calibration  tests  should  use  one  fuel  such  as  JP-4  or  JP-5,  and 
a  second  fluid  of  substantially  higher  viscosity,  for  valuation 
over  a  range  of  Reynolds  numbers  of  at  least  two  orders  of  mag¬ 
nitude.  Calibration  shall  be  conducted  near  room  temperature 
and  possibly  at  other  temperatures  such  as  50  ana  11  5  F. 

•  To  the  extent  that  scheduling  permits,  install  and  test  the  ultra¬ 
sonic  mass  flowmeter  cell  on  an  actual  full  running  engine. 
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APPENDIX  I 

SCATTERING  COMPUTATIONS 


In  this  appendix,  we  present  the  details  of  the  computations  which  were 
made  relative  to  the  acoustic  scattering  experiments.  The  scattering  par¬ 
ticles  are  in  the  quantities  given  in  Table  XIII  of  the  text  as  per  MIL-E- 
5007C.  The  objective  of  these  calculations  is  to  determine  the  total  scat¬ 
tering  cross  section  of  the  contaminant  concentration  within  a  1-cm  range 
gate  in  a  1/ 2- in. -diameter  pipe. 

We  chose  the  frequency  of  the  incident  wave  to  be  20  MHz.  A  significantly 
higher  frequency  than  this  would  have  given  rise  to  very  severe  attenua¬ 
tion  in  the  fuel,  while  a  significantly  lower  frequency  would  have  caused 
too  many  of  the  contaminant  particles  to  be  in  the  Rayleigh  region,  whose 
bound  is  given  by 

d(min)  =  ■— —  (23) 

C»  7T 

where  d  (min)  is  that  particle  diameter  for  which  geometrical  (versus 
Rayleigh)  scattering  prevails.  For  a  20  MHz  wave  and  a  sound  speed  of 
~1500  m/ s  (0.  06  in. /ps,  we  have  a  wavelength  of  12  microns  and  mini¬ 
mum  particle  diameter  given  by 

d(min)  -  12  microns  (24) 

Thus,  from  Table  XIII,  we  see  that  the  iron  oxide  particles  are  elimi¬ 
nated,  as  are  0-5  and  5-10  micron  Arizona  road  duct  particles.  The 
cotton  linters,  being  of  such  small  volume  (0.  1  gm/  1000  gal.  )  were  not 
considered.  The  particles  which  were  considered  are  listed  in  Table  XVI. 


TABLE  XVI.  CONTAMINANTS  CONTRIBUTING  TO 
SCATTERING 


No. 

Type 

Diameter 

(pm) 

Concentration 

1 

Sharp  silica  sand 

300-420 

1.  0  gm/  1000  gallons 

2 

Sharp  silica  sand 

150-300 

1,0  gm/  1000  gallons 

3 

Road  dirt 

10-20 

14%  of  8  gm/  1000  gallons 

4 

Road  dirt 

20-40 

23%  of  8  gm/  1000  gallons 

5 

Road  dirt 

40-80 

30%  of  8  gm/  1000  gallons 

6 

Road  dirt 

80-200 

9%  of  8  gm/  1000  gallons 
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The  average  diameter  of  the  particles  was  assumed  in  each  of  the  above 
instances,  e.  g.  ,  for  first  category,  we  take  the  diameter  as  360  microns. 


The  following  formula  was  derived  for  converting  the  given  concentrations 
(gm/  1000  gal.)  into  fractional  volume  concentrations: 

F  _  cone  (gm/  1000  gal,)  x  ,j_Q__6.  (25) 

v  .  ,  3.  3.785 

p  (gm/cm  ) 

3  3 

where  Fv  =  cm  of  contaminant  per  cm  of  fuel  and  p  =  mass  density  of 
contaminant. 

If  n.  is  tin  average  number  of  particles  of  the  i'  til  kind  per  cm^  of  fuel, 


i  4  3  ,3 

—  ir  a  .  it  d  . 

3  i  l 

where  d-  =  average  diameter  of  the  i1  th  particle  type.  If  d^  is  given  in 
microns,  and  n^  has  the  dimension  cm-3,  then 


n.  =  - - 

1 

it  d  . 
l 


v  12 

—  x  10 


In  terms  of  the  concentrations  listed  in  Table  XVI, 

-  =  cone  (gm/  1000  gal.)  „  2_  Jg  x  1Q5  (28) 

d  (microns  cubed) 

The  scattering  cross  section  per  particle,  assuming  roughly  spherical 
particles,  is 

0.  -  7  (d.)2  =  0.  786  (d.)2  (29) 

i  4  l 

and  for  all  particles  in  the  i1  th  category,  the  scattering  cross  section  per 
unit  volume  of  the  medium  is 


Scattering  per  unit  volume  =  n.  O’.  (• 

The  quantities  nt  and  O^  were  separately  computed  for  each  of  the  six 
particle  types.  Next,  the  overall  scattering  per  unit  volume  is  simply 
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the  linear  sum  of  the  separate  contributions: 


r  0  =  £  ni  ai  (31) 

1 

and,  finally,  if  Vc  =  volume  of  the  cell  throughout  which  we  measure  the 
scattering,  we  see  that  the  overall  scattering  is  given  by 

or  =  V  V  n.  a.  (32) 

total  c  4— '  i  i 

The  detailed  results  are  given  in  Table  XVII. 


TABLE  XVII.  RESULTS  OF  SCATTERING  COMPUTATIONS 

Particle  Type 

Average  Diameter 
(microns) 

ni 

(cm-3) 

CTi 

(cm2) 

n.  a. 
l  i 

(cm'1) 

1 

360 

4.  9  x  10”3 

1.  0  x  10~3 

5  x  10‘6 

2 

175 

4.  2  x  10‘2 

2.  4  x  10~4 

1.  0  x  10‘5 

3 

15 

74 

1.  1  x  10”6 

1.  3  x  10-4 

4 

30 

15.  5 

7.  1  x  10'6 

1.  1  x  10”4 

5 

60 

2.  5 

2.  8  x  10”5 

7.  1  x  10”5 

6 

•  ’i 

0.  44 

1.  1  x  10”4 

4.  8  x  10'5 

From  Table  XVII,  we  iir.d  that,  the  overall  scattering  per  unit  volume  of 
the  fluid  is 

^  n.  cr,  =  3.  3  x  lO-4  cm”  *  (33) 

i 

For  a  round  pipe,  the  interrogation  volume  is 

Vc  =  5  °2X  (34) 

where  D  =  pipe  diameter  and  X  =  interrogation  length.  Combining  (32), 
(33)  and  (34),  the  total  scattering  cross  section  from  the  interrogation 


90 


volume  for  a  1.  27  cm  ( 1  / 2  in.  )  diameter  pipe  and  an  interrogation  length 
of  1  cm  is 

0  ,  =  T  X  n.  a .  =  2.  86  x  3.  3  x  10  ^  £  10  ^crn^  (35) 

total  4  j  i  l 

In  order  to  check  the  sensitivity  of  the  technique  to  levels  of  this  order,  a 
glass  sphere  having  a  backscatter  cross  section  of  4  x  10"^  cm^  (diameter= 
0.  75  mm  or  0.  030  in,  )  was  lowered  into  the  fluid  (water)  and  the  Doppler 
wave  amplitude  was  measured  on  an  oscilloscope.  The  resulting  signal- 
to-noise  ratio  was  found  to  be  greater  than  30  dB.  It  was  also  found  that, 
under  flow  conditions,  when  bubbles  were  not  present  in  the  flow,  no 
Dopplers  were  observed  above  the  noise. 

Since  the  scattering  cross  section  of  a  particle  is  proportional  to  the 
square  of  its  dimension  and  the  mass  of  the  particle  is  proportional  to  the 
cube  of  its  dimension,  it  follows  from  the  above  measurements  that  the 
system  would  have  been  able  to  detect  contaminant  concentrations  of  about 
1%  of  those  given  in  MIL-E-5007C.  Such  levels,  however,  were  not 
observed  -  either  in  the  tap  water  runs  or  in  the  uncontaminated  aviation 
gas  runs.  At  this  point,  it  was  concluded  that,  even  though  the  full  two 
frequency  system  including  the  coherent  signal  processor  would  provide 
much  higher  signal-to-noise  ratio  than  that  observed  in  the  above  prelim¬ 
inary  measurements,  the  overall  accuracy  requirements  of  the  system 
are  such  as  to  raise  serious  questions  as  to  whether  the  dual  Doppler  tech¬ 
nique  is  the  correct  one  to  use  in  this  application.  The  decision,  after 
consultation  with  the  Contract  Monitor,  was  to  redirect  the  flow  velocimeter 
portion  of  the  program  in  favor  of  the  phase  measurement  described  in  the 
main  text  of  this  report.  Time  has  shown  that  this  decision  was  (at  least) 
not  unwise,  since  the  final  technique  does  in  fact  work  and  shows  high 
promise  of  leading  to  a  flow  velocimeter  which  will  accurately  measure  v 
independent  of  c. 

It  should  also  be  understood,  however,  that  while  contaminant  levels  of 
about  1%  or  more  of  those  given  in  MIL-E-5007C  appear  necessary  in 
order  for  one  to  measure  v/ c  to  1%  in  0.  5s,  with  a  given  Doppler  system, 
the  presence  of  much  lower  levels  could  be  detected. 
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APPENDIX  II 

FAR  FIELD  ANALYSIS  OF  PHASE  SHIFT  UNDER 
LAMINAR  FLOW  CONDITIONS 


PART  I.  AXIAL  WAVE 


Assume  that  we  have  a  means  of  generating  a  plane  wave  at  location  x  =  0 
in  a  pipe  of  radius  R.  The  wave  is  transmitted  upstream.  At  a  distance  x 
upstream  is  located  a  receiving  transducer, \  which  is  assumed  to  be  capa¬ 
ble  of  receiving  a  plane  wave.  Due  to  the  fact  that,  under  laminar  flow 
conditions,  the  flow  profile  is  parabolic,  rather  than  uniform,  the  wave 
traveling  upstream  will  be  transformed  from  a  plane  wave  to  one  whose 
phase  front  is  no  longer  planar.  Therefore,  the  phase  of  this  wave,  as 
measured  at  the  receiver,  will  not  be  a  linear  function  of  the  fluid  velocity. 
Simultaneous  propagation  and  reception  of  upstream  and  downstream  waves, 
and  subsequent  subtraction  of  the  two  measured  phase  shifts  will  not  cor¬ 
rect  this  problem  since  the  upstream  wave  undergoes  a  concave  distortion 
and  the  downstream  undergoes  a  convex  distortion.  Thus,  the  subtraction 
of  the  measured  phases  of  the  two  waves  results  in  0  =  (0  j  +  A)  -  ( 0 -A)  = 

(0  i  -9  2)  +  2A,  where  9  j  and  6  ^  are  the  phases  of  the  two  waves  assuming 
ideal  flow  and  A  is  the  distortion  due  to  the  actual  (parabolic)  flow  profile. 


If  we  assume  that  the  wavelength,  X  0,  of  the  acoustic  wave  is  much  smaller 
than  the  pipe  radius  R,  and  that  the  distance  x  between  transmitter  and  re¬ 
ceiver  is  large  enough  so  that  x  >>  R^/X  D,  we  can  rather  easily  obtain  the 
expression  for  the  measured  phase  at  point  x.  We  first  divide  the  cross- 
sectional  area  of  the  pipe  into  a  number  N  of  concentric  rings,  each  with  a 
width  of  AR.  The  measured  phase  will  be 


0 


N 


V.  2tr  R,  AR  sin  0 
.  -1  ;  i=i  1 

tan  \  — jq - 

Y'  2rr  R.  AR  cos  0  . 

k  1 


A 

r| 


(36) 


where  R^  =  radius  of  i'th  strip  and  0  ^  =  phase  of  the  i'th  wave  component  at 
point  x: 


0 


Zt  x 

X  . 

1 


(37) 


where  X  ^  =  wavelength  at  the  radius  R^  . 


Equation  (36)  can  be  written  in  integral  form  by  letting  AR^-dr: 
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R 


<J>  =  tan 


..rX 


r  sin  0  (r)  dr 


~F~ 

i tt  /  r  c< 


2 it  /  r  cos  0  (r)  dr 
'0 


whe  re 


0  (r)  = 


2  TT  X 

^  (r) 


If  the  flow  is  laminar,  the  flow  profile  will  be  parabolic,  i.  e.  , 

v(r)  =  v(0)  [  1  -  (“)2] 
where  v(0)  “  flow  velocity  at  r  -  0. 

Therefore,  since 

Mr)  = 


where  f  transmitter  frequency, 


Mr)  = 

from  which  we  find  0  (r): 

0  (r)  = 


c  -  v(0)  f  1  -  (£f] 


_ ux _ 

c  -  v(0)  [  1  -  (■£)  ] 


The  overall  phase  of  the  received  wave  is,  from  (38), 

R 


0  =-•  tan 


l  ?  .  r  wx  1 

I  /  r  sin - 

-l)J0  [c  -  v(0)  [  1  -  (|)  ]_ 


r  r  wx 

f  r  cos - - - 2~ 

0  Lc  -  v(0)  [!-(“)]. 


dr 


dr 


0  c  -  v(0)  [ 

Note  that,  when  v(0)  =  0,  the  above  equation  has  the  solution 

0  (v(0)  =  0)  =  ~ 
c 

which  is,  of  course,  the  correct  solution  for  zero  flow. 


Equation  (44)  can  be  solved  by  transforming  the  variable  from  r  to 
given  by  Eq.  (43). 


(38) 


(39) 


(40) 


(41) 


(42) 

(43) 


(44) 


(45) 

0  (r)  as 
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Abbreviatinu  Rq.  (44)  am 


0  «  tan 


"It! 


wo  find  that 


ami 


tl2 

R  wx 

il  <R) 

r 

_L_ 

a  | |i  /I 

A  il 

2  v(Q)  J 

II  (0) 

mu  f* 

U  fr 

ni 

R  wx 
*  2  v(0) 

7 

rt  (0) 

o 2 

cott  0  d  0 

the  solutions: 

R2  ux  { 

1  sin  • l 

an 

.  f 

COS  l> 

ddj 

0  (R) 

2  v(0)  j 

1  0 

1) 

0 

0  (o) 

R^  ux  j 

( cos  0 

1) 

.  f 

sin  (1 

d(>  j 

0  (R) 

2  v(0)  | 

0 

0 

n  (o) 

Thus,  the  measured  phase  at  point  x  ulonn  the  pipe  is 

0  <K) 


0  --  tan 


-1 


sin  (>  ( R )  sin  0  (0)  f  cos  0 
~0  (R)  "  (MO)  0 


jam. 


cos  0  (R)  cos  0  (Qj 

(MR)  "  «(0) 


0(0) 

L'J  it1 


0  (R) 


(4(») 


(47) 


(4HI 


(40) 


(50) 


(51) 


Wo  next  note  from  Eq.  (43)  that  0  (R)  is  just  the  phase  under  aero  flow  con¬ 
ditions.  Thorofore,  0  (R)  can  bo  fixed  and  will  romain  constant  for  any 
(laminar)  flow  conditions.  For  convenience,  wo  set 


0  (R)  a  2nir  ,  n  ■-  1,2,3 . 

This  means  that  0=0  under  zero  flow  conditions, 

We  now  have  a  simpler,  but  no  less  exact,  solution  for  0  : 


(52) 
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0(R) 

Aituilfll  .  f  mJ 

0  (0)  J  (1  H 

.  .  ■  .1M _ — _ 

0  (0) 

zz&JLM  .  -J.-  .  /  sinJ 

0(0)  0(R)  /  « 

0  (R) 

We  next  wish  to  determine  whether  nr  not  Rq,  (51)  behave*  linearly  in  the 
relationship  between  0  and  v(U).  Although,  strictly  speaking,  numerical 
method*  ahould  he  employed.  If  we  restrict  our  solution  to  small  value*  of 
0  (0)  •  0  (HI,  we  can  obtain  a  measure  of  the  range  of  linearity  of  0  v*  v(Q), 

We  first  note  that  0  (0)  i*  large  because  the  total  path  contain*  many  wave¬ 
lengths,  But  if  we  restrict  the  range  0  (0)  -  0  (R)  so  that 

|  0  (0)  -  0  (R)!  <  <  ~  ,  we  can  make  the  following  approximation*) 


sin  0  (0) 

*  0  (0)  . 

■  0  (R) 

(64) 

cos 0  (0) 

r.  J  .  1 

2 

1 0  (0)  -  0  (R))  2 

(55) 

a 


l  -  |  [  0  (0)  -  0  (R)2  -  1 
\  [0  (0)  4  0  (R)] 


[0  (0)  -  0  (R)] 
[  0  (0)  +  0  (R)J 


(56) 

(57) 

(58) 


The  justification  for  (56)  is  that  0  »  2ntr  over  the  integration  range.  There¬ 
fore,  cos  0  =  1  over  this  range.  The  opposite  state  of  affaira  occurs  with 
(57),  whero  we  notice  that,  since  0  -  2nrr  is  large  (as  we  originally 
assumed),  and  since  the  integration  range  is  small,  we  can  allow  sin  0  to 
be  the  variable,  and  in  the  denominator  let  0— »(l/2)  [0  (0)  +  0  (R)]  . 


Therefore,  under  the  conditions  specified, 
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0 


(  ln(  ) 
un-1  )  — - HP!..  l(M0)/ 


'  0  (0) 


(0)  •  0  (R)| 


-J_  *  l.tf Iftl  -lMHH 

f‘  I  I)  (0)  t  0  (H) I 


I 


Hiiu-'o 


1  TTto)  -lui  -HToTMMHi) 

Kcji  (M)  becurneti 

( !!  19)  jHJ  .  ln  (£1111)  ) 

0  ;  un-*l _ *M _ -  V<M0)  I  ( 

)  >  1 

(  0(G)  *  (MR)  ) 


0  (0)  (MR) 
lotting  p'(0)  -  v(0)/c,  aiul  noting  that 


aiul 


0(0)  ••= 


(MR)  ; 


ux 


c  -  v(0) 


w\ 


wo  fine!  that 


0  (0)  -  0  (R)  -  p  '(0) 


ln  (fyo))  ;  ln,ltP' (0) )  ’  P  *  (o)  -  i  P 


(0) 


1 


l 


0  (0)  0  <R)  wx 

Substituting  the  above  three  quantities  into  (61),  we  finally  obtain 

-1  )  n  tr_x  /v(0)\ 


.  .  *1  nirx  „ 

<P  *  tan  /  ~ —  (1  '(0) 


tan 


(*‘0 


0>Q) 


(61) 


(6*) 

(63) 


(64) 

(65) 


(66) 


For  parabolic  flow  profile  the  average  velocity,  v  is  just  half  the  maximum 
volocity  v(0),  Therefore, 


<P  =  tan 


- 1  j  2jvit x  /  v  \  I 


(67) 


96 


The  approximations  required  to  obtain  this  relationship  worn  accurate 
to  order  p"*,  Therefore,  wo  find  that  (f»?)  la  linear  in  v(G)  only  inaofar 
aa  the  arc  tangent  la  linear.  For  example,  If  inaccuracy  of  J_  0,  t,$%  were 
require*!,  the  ranye  of  0  could  not  exceed  4.  8  dey  and  the  meaavirement 
preciaion  of  0  itaelf  would  have  to  he  better  than  £  0,  12  dey!  If  inaccuracy 
of  +_  0,  6%  were  required,  the  allowed  ranye  of  0  would  bo  9,  i  .ley  and  a 
meaaurement  preciaion  on  0  would  be  ±  0.  Ob  dey, 

To  extrapolate  the  interpretation  of  (u7)  further  would  be  to  violate  the 
aaaumptiona  made  in  deriviny  it.  However,  it  la  obvioua  that  Rq,  ( 5  i ) 
(which  ia  the  exact  aolution  for  the  problem  formulated)  becomea  inereas- 
inyly  nonlinear  aa  the  ranye  0  (0)  -  0  (R)  ia  widened.  It  seems  reaaonable 
to  conclude  that,  within  the  potential  ranye  of  hiyh  accuracy,  the  instru- 
mental  errora  in  mcaaurlny  phase  angle  (to  within  ~  0,  06  deg)  would 

prevent  high  accuracy  from  actually  being  achieved, 

It  la  clear  that  instrumental  errors  prohibit  high  accuracy  for  small  values 
of  0  ,  and  the  nature  of  the  flow  itaelf  prohibits  high  accuracy  for  large 
values  of  0  , 

It  ia  precisely  because  of  these  limitations  on  the  usual  "upstream  minus 
downstream"  phase  or  time  difference  flowmeter,  that  we  considered  a 
dual-frequency  measurement  based  on  scattering. 


PART  H.  OUHQUE  WAVE 


We  wish  to  determine  if,  to  first  order  approximation,  the  measurement 
of  phaso  is  independent  of  flow  profile  in  an  oblique  transmission  systom. 
We  assume  piano  wave  transmission  and  reception.  The  sotup  is  shown 
in  Figure  41.  A  plane  wave  is  transmitted  to  a  plane  wave  receiver  across 
a  circular  pipe  of  radius  R.  The  beam  makes  an  angle,  p  ,  with  respect 
to  the  normal  to  the  axis  ol  the  pipe.  If  the  receiver  and  transmitter  trans¬ 
ducer  apertures  are  very  small  as  compared  with  the  radius  R,  then  the 
phase  0  of  the  received  wave  is  giver,  by 


where 


0  = 

f  2tt  dL 

J0  X. 

(68) 

1 

2  R 

(69) 

Li 

CU3  P 

clL  = 

2  dr 

(70) 

cos  P 
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Figure  41.  Wavefront  Distortions  for  Oblique  Incidence 
Laminar  Flow. 
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In  the  ertNti  of  a  parabolic  flow  profile, 


and 


X.  .  Mr)  *  2  ”  •  ''!>» 

y 

v(r)  -  v(0)  f  1  “  (  R  >2  ) 


(71) 

(72) 


where 


v(r)  "  flow  velocity,  v(0)  =  axial  flow  velocity, 
c  =  sound  speed,  and  w  =  angular  frequency. 


Substituting  (69),  (70),  and  (71)  into  (68)  yields 

R 
to 

<P 


?.  to  r 

c  cos  P  JQ  J 


dr 


P’(0)[  1  -  (“)2] 


(73) 


where 


P  MO)  ev(0)/c 


(74) 


Assuming  that  p'(0)  <<  1, 


R 

*  i  r^rp  +  -  <£>2 1}-< 


(75) 


The  change  in  phase  due  to  a  change  from  zero  flow  to  finite  flow  is 

R 


*  i  f  [i.  (|)2]  a 

r  C  COS  P  Jq  R  1 


(76) 


The  total  mass  flow  rate  in  the  pipe,  M,  is  given  by 

R 


M  =  f 


M  =  /  p  v(r)  dA 
'0 

where  p  =  mass  density,  assumed  constant,  not  a  function  of  r 

From  Eq.  (72)  and  noting  that  dA  =  2  tt  rdrwe  have 

R 


L 


(77) 


ivt  =  2ir  p  v(0)  r  [  1  -  (-)  ]  dr 


(78) 


or 


R 

Kl  =  2TrpcP'(0)  /  r  [  1  -  (~  )  J  dr 

Jq  R 
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Thun,  the  ratio  of  phase  shift  to  flow  rate  is 


Jvi 


R 

f  ,  t .  (x  j2 ,  dr 

c  cos  3  Jo  1  1  R  J  tlr 


2  np  c  p  '(0) 


R 

l 


r  [  1  (^)2  ]  dr 


(79) 


Wo  soo  in  the  above  equation  that  the  axial  flow  term  p'(0)  cancels,  and 
that  AO/M  is  thus  independent  of  (laminar)  flow  velocity  profile  v(r).  That 
is  to  say,  the  measured  phase  shift,  or  time  interval,  is  proportional  to 
mass  flow  rate  Ivt,  for  laminar  flow,  when  the  stated  assumptions  am  met. 
Instrumental  calibration  is  required,  however,  for  different  pipe  sizes 
because  A</)/h/l  is  obviously  dependent  upon  R. 


It  should  be  understood  that  although  the  wavefront  is  assumed  plane  at  the 
transmitter  and  receiver,  the  wavefront  becomes  curved  and  then  uncurved 
as  it  traverses  the  fluid,  as  far  as  its  projection  in  the  plane  of  incidence 
is  concerned.  Actually,  however,  its  projection  increasingly  curves,  in 
planes  orthogonal  to  the  plane  of  incidence,  and  parallel  to  the  axis. 


The  above  remarks  apply  to  the  far  field.  In  the  near  field,  to  the  extent 
that  waves  remain  plane,  phase  measurements  are  essentially  proportional 
to  the  How  velocity.  Referring  to  Figure  41  again,  we  may  note  that  if  one 
dimension  of  the  square  flow  channel  is  modified,  giving  a  rectangular 
cross  section  of  different  area,  one  can  arrange  to  extend  or  narrow  the 
flow  range,  for  a  given  set  of  fixed  electrical  parameters  such  as  frequency 
and  range  of  phase  shift. 
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